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Abstract. We study the large-A'' limit of a system of bosons in- 
teracting with a potential of intensity 1/A'^. When the ground state 
energy is to the first order given by Hartree's theory, we study the next 
order, predicted by Bogoliubov's theory. We show the convergence of 
the lower eigenvalues and eigenfunctions towards that of the Bogoliubov 
Hamiltonian (up to a convenient unitary transform). We also prove the 
convergence of the free energy when the system is sufficiently trapped. 
Our results are valid in an abstract setting, our main assumptions being 
that the Hartree ground state is unique and non-degenerate, and that 
there is complete Bose-Einstein condensation on this state. Using our 
method we then treat two applications: atoms with "bosonic" electrons 
on one hand, and trapped 2D and 3D Coulomb gases on the other hand. 
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1. Introduction 



In a famous paper |9j, Bogoliubov was able to predict the excitation spec- 
trum of a quantum gas satisfying the Bose statistics and he used this to 
understand its superfluid behavior. Since Bogohubov's work, there has 
been several attempts to formulate Bogoliubov's theory in a mathemati- 
cally rigorous way. This was especially successful for completely integrable 
ID systems [211 [351 [32l [12l [TTl [591 [60] , for the ground state energy of one and 
two-component Bose gases |41l 1^21 [55] , and for the Lee-Huang-Yang formula 
of dilute gases [191 1221 [H] • Recently, Seiringer [52] and Grech-Seiringer |24] 
have for the first time justified Bogoliubov's theory for the excitation spec- 
trum of trapped Bose gases, with a general short range interaction, in the 
mean-field regime. See, e.g., [65| for a recent review on the subject and [15] 
for a discussion of translation- invariant systems. 

The purpose of this article is to give general conditions under which Bo- 
goliubov's theory is valid, that is, predicts the lowest part of the spectrum 
of the many-body Hamiltonian of bosons, in the mean-field regime. Our 
results cover a very large class of interacting Boses gases and they general- 
ize the recent works |52[ I24|. In particular, our method applies to Coulomb 
systems. 

We consider a system of N quantum particles, described by the Hamil- 
tonian 



of square-integrable functions ^' G L^(il^) which are symmetric with respect 
to exchanges of their variables, namely 



for every a in the permutation group ©at. Here Q is an open subset of 
with d > 1, T is a self-adjoint operator on L^(r2) with domain D{T), 
and w is an even real-valued function describing the interactions between 
particles. We have neglected spin for convenience, but it can be added 
without changing any of our result. The Hamiltonian describes a system 
of bosons living in il.. 

The operator T can contain both the kinetic energy and an external po- 
tential which is applied to the system, including possibly a magnetic field. 
We typically think of T = —A on a bounded set Q with appropriate bound- 
ary conditions (Dirichlet, Neumann or periodic), or of T = —A + V{x) on 
Q = R*^, with V an external potential which serves to bind the particles. 
In the latter case, the function V could tend to zero at infinity but it then 
has to be sufficiently negative somewhere, or it could tend to infinity at 
infinity, in which case all the particles are confined. We could also replace 
the non-relativistic operator —A by its relativistic counterpart \/l — A — 1. 
We shall keep the operator T sufficiently general in this paper, such that all 




i=l ^<i<j<N 

acting on the symmetric (a.k.a. bosonic) space 

N 



iD^ = (g)L2(J7) 



sym 
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these situations are covered. The function w could also be replaced by an 
abstract two-body operator but we do not consider this here for simplicity. 

We are interested in the limit of a large number N of particles. Here we 
are considering the mean-field regime, in which the interaction has a fixed 
range (the function w is fixed) but its intensity is assumed to tend to zero 
in the limit — )• oo, hence the factor 1/{N — 1) in front of the interaction 
term in the Hamiltonian H^- This factor makes the two sums of order N 
in Hn and, in this case, an important insight is given by Hartree theory. 

Let us recall that a Hartree state is an uncorrelated many-body wave 
function in which all of the particles live in the same state u G L^(il) such 
that |np = 1, and which takes the form 

^'(xi, ...,xn) = u{xi) ■ ■ ■ u{xn)- 

The energy of such a state is 

{^,Hn^) = N ({u,Tu) + ^D{\u\^,\u\^)] :=N£n{u) 



where 



Dif,g):= / / f{x)g{y)w{x-y)dxdy 
Jn Jn 



is the classical interaction. Henceforth, all the Hilbert spaces we consider 
have inner products which are conjugate linear in the first variable and linear 
in the second. 

Provided that there is Bose-Einstein condensation, the leading term of 
the ground state energy 

E{N) := inf spec H]\f 
is given by Hartree's theory: 



E{N) = Nen + oiN), 



where en is the corresponding Hartree ground state energy: 

en:= inf £niu) = inf \ {u,Tu) + lDi\u\Mu\^)] . (1) 

||n||=l |H|=1 

In this paper, we shall assume that there exists a unique Hartree minimizer 
uq for en. It is then a solution of the nonlinear Hartree equation 

= {T +\uo\'^ *w - tin)uo := huo, (2) 

where //h G K is a Lagrange multiplier. 

Bogoliubov's theory predicts the next order term (of order 0(1)) in the 
expansion of the ground state energy E{N). It also predicts the leading 
term and the second term for the lower eigenvalues of H^- The Bogoliubov 
method consists in describing variations of the wavefunctions around the 
Hartree state uq ■ ■ ■ uq in a suitable manner. We will explain this in 
detail in Section 12.21 below. The final result is an effective Hamiltonian H, 
called the Bogoliubov Hamiltonian which is such that the lower spectrum of 
Hj\[ in S)'^ is given, in the limit — )• oo, by the spectrum of the effective 
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operator Neu + IHI, up to an error of order o(l). Formally, we therefore find 
that 



Hn ^ Neu + B[ + o{l). (3) 



This vague statement is made precise in our main result, Theorem [2] below. 

The essential fact about the Bogoliubov Hamiltonian HI is that it is a non 
particle conserving self-adjoint operator acting on the Fock space 

n 

^+:=C©0(g)i^+ 

n>l sym 

where 

^+ = {uo}^ c L\n) 

is the one-body space of excited particles. We started with a particle- 
conserving model and we end up with a theory in Fock space, in which 
the number of particles is not fixed. The reason is that we are describing 
here the excitations around the reference Hartree state. 

For the acquainted reader, we mention that HI is indeed nothing but the 
second-quantization of (half) the Hessian of the Hartree energy at uq, and 
its expression in a second quantized form is 



/ a*{x){ha){x)dx +^ [ [ w{x - y)(2uo{ 
Jn ^ JnJQ \ 



x)uo{y)a*{x)a{y) 

+ uo{x)uo{y)a* {x)a* (y) + uo{x)uo{y)a{x)a{y)^ dx dy. (4) 

Here a*{x) is the creation operator of an excited particle at x, acting in the 
Fock space J-+ , and h is defined in ^ . We will explain the meaning of this 
formula later in Section 12.21 

A result similar to ([3]) has recently been obtained for weakly interacting 
Bose gases by Seiringer [52] and Grech-Seiringer [23]. They assumed that 
w is bounded, decays fast enough and has non-negative Fourier transform. 
The operator T was T = —A in a box with periodic boundary conditions 
in [52] and T = - A + V{x) on R'^ with V{x) +oo at infinity in [23]. Our 
method is different from that of |52l and it applies to a larger class of 
models. 

We give in this paper a list of abstract conditions that a Bose gas should 
satisfy in order to get the Bogoliubov result These conditions are given 
and explained in Section [2. II below. Loosely speaking, we assume that there 
is complete Bose-Einstein condensation on a unique Hartree minimizer uq 
which we assume to be non-degenerate. Our message is that, once the Bose- 
Einstein condensation is proved, one can get the next order in the expansion 
of the energy by Bogoliubov's theory. No further assumption is needed. 

The paper is organized as follows. In the next section, we define our model 
by giving the appropriate assumptions on T and w, and we properly define 
the Bogoliubov Hamiltonian M. We then state our main results, Theorem [2] 
and Theorem [3l In Sections 13.11 and [3T2l we apply our abstract result to two 
particular examples: bosonic atoms and trapped Coulomb gases. Sections[5]- 
[7] are devoted to the proof of the main abstract results. 
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2. Main abstract results 
In this section we state our main result. 

2.1. Assumptions. We start by giving the main assumptions on T and w, 
under which our results apply. Later in Section [3] we consider two specific 
examples, for which all the following assumptions are satisfied. 

The first condition concerns the properties of T and w which are necessary 
to give a proper meaning to the many-body Hamiltonian Hi\i. 

(Al) (One- and two-body operators). The operator T : D{T) — )- L'^{Q) is 
a densely defined, bounded from below, self-adjoint operator. The function 
If : M'^ — )• M is Borel-measurable and w{x) = w{—x). Moreover, there exist 
constants C > 0, 1 > ai > and a2 > such that 

- ai(T, + Ty + C)< w{x -y)< a2(r, +Ty + C) on L^{n''). (5) 

Note that, although we keep the one-body operator T abstract, we use a 
two-body operator w which is a translation-invariant multiplication operator 
in L^(f]^). This is only for convenience. All our results are also valid if 
w is an abstract two-body operator on S)"^ which satisfies an estimate of 
the same type as ([5]), and if = L?'{Q) is an abstract separable Hilbert 
space. However, in this case the expressions of the Hartree energy and of the 
corresponding nonlinear equations are different (they cannot be expressed 
using a convolution). We shall not consider this abstract setting to avoid 
any confusion. 

Under Assumption (Al), is bounded from below, 

N 

HN>{l-ai)Y,T,-CN. (6) 

i=l 

In the paper we always work with the Friedrichs extension [l6] , still denoted 
by Hn. Note that we do not assume the positivity or boundedness of w 
or its Fourier transform, but only that it is relatively form-bounded with 
respect to T. 

Our second assumption is about Hartree theory. 

(A2) (Hartree theory). The variational problem (d) has a unique (up to a 
phase) minimizer uq in the quadratic form domain Q{T) of T . Moreover, 
uq is non- degenerate in the sense that 

[ K* h + Kl J ^^^^ ^+®^+ 
for some constant r]ii > 0, where io+ = {no}"*". Here 

h := T + |uo|^ * u) — //H) 
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with the Lagrange multiplier fin := eii + D{\uo\'^, \uq\'^)/2 (which ensures that 
Huq = 0), and Ki : S)j^ — )• i3+ and K2 = — )• io+ are operators defined by 



{u,Kiv) = u{x)v{y)uo{x)uo{y)w{x - y)dxdy, 

Jn Jn 



{u,K2v) = u{x)v{y)uo{x)uo{y)w{x - y)dxdy 

Jn Jn 

for all u,v £ Sjjf- . The operators Ki and K2 are assumed to he Hilhert- 
Schmidt, that is 

\uo{x)\'^\uo{y)\'^w{x - yfdxdy < 00. (8) 

n Jn 

Remark 1. Note that once we have assumed that uq is a minimizer for 
then we always have the mean-field equation Huq = due to a standard 
argument (see e.g. [36, Theorem 11.5]j. In the whole paper we will for 
simplicity use the same notation for the operator h on the full one-body 
space Sj and for its restriction to the smaller space Sj^. For an operator 
A on Sj^, the notation A means JAJ with J being the complex conjugate, 
namely A{v) = A(v) for every v G S)j^. 

Remark 2. While we shall treat Ki as a one-body operator, we should really 
think of K2 as its integral kernel K2{x,y) = [Q Q) {uq uow{. — .)){x,y), 
which is the two-body function obtained by projecting the symmetric function 
Uoix)uQ{y)w{x - y) onto S)\. 

In (A2) we are making assumptions about the uniqueness and non-degeneracy 
of the Hartree ground state mq- The Hessian of the Hartree energy can easily 
be seen to be 

^Hess £n{uo){v,v) 
1 



2 



n Jn 



{v,hv) + - / w{x-y)(v{x)uo{x)uo{y)v{y)+v{x)uo{x)uo{y)v{y) 



+ v{x)uo{x)uo{y)v{y) + v{x)uo{x)uo{y)v{y) j dx dy 

= 2((S)'CKr h+ir^(l))^^^-^ (9) 

for all V G It turns out that the non-degeneracy of the Hessian, 

Ress £uiuo){v,v) > 2r/H ||f 11^2(0) , 

is equivalent to our assumption ([7]). When uq is real, as it is in many 
applications, then io+ = S)^, Ki = K2 and it can be verified (using a test 
function of the form —v)) that ([7]) implies h > rjn on S)^, which means 
that there is a gap above the first eigenvalue 0. In general, we however only 
know that h -\- Ki > t^h • 

The Hilbert-Schmidt assumption ([8|) on Ki and K2 will be useful later to 
ensure that the Bogoliubov Hamiltonian H is well defined (see Section | 
below). 
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Our last assumption is about the validity of Hartree theory in the limit 
— )• oo. We assume that the system condensates in the unique Hartree 

ground state uq. This assumption will be necessary for the proof of the 

lower bound on the spectrum of H^. 

(A3) (Complete Bose-Einstein condensation). For any constant R > 0, there 
exists a function ej? : N — )• [0, oo) with liniAr^oo £_r(-^) = such that, for 
any wave function G Sj^ satisfying {"^NjHn'^n) < E{N) + R, one has 

i^^°^^>l-.MA^) (10) 

where uq is the Hartree minimizer in Assumption (A2). 

Here 7^ is the one-body density matrix of the wave function ^' E S)^ , 
which is the trace-class operator on L^(r2) with kernel 



7^(x,y):=iV J 'i/{x,X2, .■.,XN)'^{y,X2, ■.■,XN)dx2...dxN- 



Note that a Hartree state has the density matrix j^^n = Nu{x)u{y). There- 
fore ([Toll is the same as saying that 741 is in some sense close to 7„8iAf- For 
more explanation about the Bose-Einstein condensation, we refer to the 
discussion in [39]. 

In many applications (especially for Coulomb systems, see Section [3]), a 
stronger condensation property will hold true. Namely, we will have a bound 
from below valid for all G i^^, and not only for those which have a low 
energy. We therefore introduce the following stronger assumption, which 
obviously implies (A3): 

(A3s) (Strong condensation). We have h > r]}i > on Sj^, and there exists 
a constant < eg < 1 such that 

N 

Hn - Nen > (1 - eo) Yl + °(^)- 
i=i 

Here h is the mean-field operator given in Assumption (A2). 
In fact, in practice (A3s) follows from a Lieb-Oxford inequality 

/^A Wij]^) > Id{p<!,,p^) + error 



where pqi^x) = 7iir(x,x). It is for proving estimates of this form that it is 
often useful to know that w > where ^ denotes the Fourier transform. We 
will come back to this in Section [3] where we consider two examples. 

2.2. The Bogoliubov Hamiltonian. Near the Hartree minimizer we have 
■^^( (1 + 11^12)1/2 ) = £u{uo) + \iiess £u{uo){v,v) + o{{v,iT + C)v)) (11) 

for any v which is orthogonal to uq, that is, v £ The next order in 

the expansion of the eigenvalues of the Hamiltonian Hj^ will be given by 
the Bogoliubov Hamiltonian H, which is obtained by second quantizing the 
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Hessian in (jlip . More precisely, this means replacing v{x) by an operator 
a*{x) which creates an excited particle at x, and v{x) by an operator a{x) 
which annihilates it. These operators (formally) act on the Fock space of 
excited particles 

n 

^+ :=ce0(g)io+ = ce0i3:^. 

n>l sym n>l 

So the expression of the Bogoliubov Hamiltonian is 
a*{x){{h + Ki)a){x)dx 

1 



2 



Q. Jn 



K2{x,y)a*{x)a*{y) + K2{x,y)a{x)a{y) ]dxdy. (12) 



In order to make the formula (I12p more transparent, let us explain how 
the Hamiltonian EI acts on functions of J^-^ . If we have a ■i/'fc £ -5+, with 
k > 2, then we get 

HV'fc = • • • e © e j ^{h + Ki)j I Vfc e e VI.+2 © • • • (i3) 

where 

il:'k+2{xi,...,Xk+2) =^=^==== ^ ^2(2;<T(i),a:<^(2))V'fc(a:<x{3),---,2;^(fc+2)), 



ilj'k_2{xi,...,Xk-2) = \/k{k-^) I dxfc_i / dxfcK2(a;fc-i,2;fc) V'fe(a;i,---,a;fc)- 

JO. 

The link between the formal expression ()12p and the rigorous formula ()13p 
is explained in |3T1 Sec. 1]. See also ([67|) in Appendix |A] for another equiv- 
alent expression of HI using one-body density matrices. 

Let us remark that for ™ L^(r2'^+^) for ah G L'^{VL^), it 

is necessary and sufficient to have K2{.i ) in L?'{VL^), which is the same as 
assuming that K2 is a Hilbert- Schmidt operator, as required in Assumption 
(A2). 

Since Ki and K2 are Hilbert- Schmidt, the Hamiltonian EI is well defined 
on states living in truncated Fock spaces and in the domain of h: 

M n 

U 0(2)^w- (14) 

M>0 n=0 sym 

The following theorem tells us that EI is bounded from below and it is a 
well-defined self-adjoint operator by the Friedrichs method with the form 
domain being the same as that of dr(l + h) on Here we have used 
the usual notation dr(A) for the second quantization in Fock space of an 
operator A acting on the one-body space S)+: 

00 m „ 

dr(^):=0^^j= / a*{x){Aa){x)dx. (15) 

m=0 i=l 
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Remark 3. We shall always denote by C > some (large) constant which 
depends only on T and w . Two C 's in the same line may refer to different 
constants. 

Theorem 1 (Bogoliubov Hamiltonian) . // (A1)-(A2) hold true, then the 
operator M is symmetric and on the core domain ()14p one has 

C-^dT{h + 1) - C7 < M < dr{h + C) + C. (16) 

Consequently, the form domain of its Friedrichs extension (still denoted by 
M) is the same as that o/dr(l + h) on T^: 

oo n 

0(g)I)((/i + l)V2). (17) 

n=0 sym 

Moreover, we have the following spectral properties. 

(i) (Ground state and ground state energy). The Hamiltonian M has a 
unique ground state in Tj^. The ground state energy o/H is strictly negative, 
except ifw = 0, in which case we have mia{M) = (the ground state being 
the vacuum in ). 

{a) (Essential spectrum). The essential spectra of h and H are simultane- 
ously empty or non empty, and we have in the latter case 

aess(IH) =a(]H)+aess(/i). 

Consequently, inf cJcss(H) — inf cr(]HI) = inf (Jcss(/i) > ??H > 0. 

[Hi) (Lower spectrum). Assume that T = T (in this case uq is a real-valued 
function and hence Ki = K2). If h + Ki has infinitely many eigenvalues 
below its essential spectrum, then H also has infinitely many eigenvalues 
below its essential spectrum. 

On the other hand, ifT = T, Ki > and h has only finitely many eigen- 
values below its essential spectrum, then H also has finitely many eigenvalues 
below its essential spectrum. 

We refer to Appendix lAl for a proof of Theorem [1] and further discussions. 

Remark 4. Note that since h+Ki > rjn due to ^ andKi is Hilbert- Schmidt, 
we always have inf acssih) > f?H > 0. 

Remark 5. The reader should be cautious with the fact that, when w ^ 0, 
even though K2 is a Hilbert- Schmidt operator on the one-body Hilbert space, 
the pairing term 

^ Inin (^2(a;,y)a*(2;)a*(y) + K2{x,y)a{x)a{y)^ dx dy 

is neither bounded on J^^, or relatively compact with respect to dT{h). In- 
deed, when the essential spectrum of h is non empty, we have 

inf fTess(IHI) = inf fT(EI) +miacss{h) < inf a ess (h) = inf cTess(dr(/i)) 

due to TheoremUl and hence H and dT{h) have different essential spectra. 
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2.3. Exciting the Hartree state. In this section we explain how to ap- 
propriately describe the variations of a many-body wave function ^ in the 
neighborhood of the Hartree state Uq'^ . We will see that the Fock space 
arises naturally. We will also define a unitary operator Un which will be 
essential for our main statements. 

For two symmetric functions *fc G Sj'' and G Sj^, we recall that the 
symmetric tensor product is defined by 

Note that this tensor product satisfies the commutative property ^i^-f^s^e = 
^i^s'^k and the associative property {'i>k0s^e)®s'^m = *ik®s (*£®s^m) 
for all ^ ^ and £ Consider now any (real- valued) 

orthonormal basis uq, wi, ... of containing the Hartree minimizcr uq. 

Then, it is known that {ui^ ®s ' " ®s ^ijvl is an orthogonal basis of the 
symmetric space , where 



ui^s---®sUn{xi,...,xn) ■■= -y^ ^ u^(^i){xi)---u^(^n){xn)- 



o-g6jv 



So we can write our many-body Hilbert space as a direct sum of spaces 
where ^ = span(?xo (8) • • • (8) uq) and 



k 

N-k 

where we recall again that Sj+ = {uq}^ = spanjui, ?i2, •••} C L'^{Cl). In 
other words, we can write any wavefunction ^ G S)^ as follows 

* := ^0 ur + uT~'^ ®s ^1 + nf''-'^ ^2 + • • • + ^iv 
where ipk € S}^. It is a simple exercise to verify that 

from which we deduce that 

N 

k=l 

Therefore we see that the linear map 

JV 

C/iv: ^ ^|^ = 0i5!; 

n=0 

* 1-^ V'O © V'l ® • • • © '0Ar 

is a unitary operator from S)^ onto the truncated Fock space J='^'\ The 
latter can always be seen as being embedded in the full Fock space T'^ of 
excited particles and the unitary operator C/jv is also a partial isometry from 



(18) 
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i^^ to . We see that, in the Umit N ^ oo, the Fock space of excited 
particles arises naturally as the hmit of the truncated Fock spaces J-^^ . 

The operator Un is a mathematical tool which implements what is called 
a c-number substitution [401 139j. In Fock space the usual way to formalize 
the c-number substitution is to use the Weyl operator, and here Un plays 
the same role. The difference is that the Weyl operator is defined on the 
whole Fock space J- with values in J^, while with the operator C/^v we go 
immediately from the A/^-body space S)^ to the excitation Fock space 
which is a proper subspace of T. We shall give important properties of the 
operator Un in Section HI 

One of our main results will be that 

in an appropriate (weak) sense. 

2.4. Convergence of the excitation spectrum. A convenient way to 
describe the lower eigenvalues is to use the min-max principle (see [37]). If 
^ is a self-adjoint operator, which is bounded from below, on a (separable) 
Hilbert space, then we may define the min-max values 

Ai(^) := inf max($,^$) 

Y subspace <E>g5^ 
dimY=L \W\=1 

for L = 1,2, .... It is known that limi_>>oo ^li^^) = inf (Jess ("4.), where we use 
the convention that inf cress(-^) = +oo when the essential spectrum of A is 
empty. Moreover, if 

Xl{A) < inf CJess(^), 

then {Xj{A)}j^i are the lowest L eigenvalues of ^, counted with multiplicity. 
Our main result is the following. 

Theorem 2 (Convergence of the excitation spectrum). Assume that (Al)- 
(A2)-(A3) hold true. 

{i) (Weak convergence to H) . For every fixed $ and in the quadratic form 
domain of the Bogoliubov Hamiltonian H, we have 



JV— S.OO + + 



(19) 



where Un is defined in and by convention is extended to outside 
ofT^ . 

(ii) (Convergence of eigenvalues). Let Xi{H]\[) < X2{H]\[) < ... and Ai(IHI) < 
A2(]HI) < ... be the min-max values of in and IH in respectively. 
We have 



lim [XL{HN)-Nen] = Al(H) 

N^oo 



for every L = 1,2,.... Consequently, we have the spectral gap 
liminff inf (Jess (^jv) - Xi{Hn)) > inf aess(IH) - Ai(EI) = mfa^^h) > 0. 

N->-oo \ J 

(Hi) (Convergence of the ground state). The Hamiltonian H has a unique 
ground state ^^^^ in (up to a phase factor) and also has a unique 
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ground state ^'j^^ for N large enough. Furthermore (up to a correct choice 
of phase for ^-^^J, 



lim Un'^m^ = ^^^^ 



(20) 



strongly in The latter convergence is strong in the norm induced by the 
quadratic form of M. on F+ if (A3s) holds true. 

{iv) (Convergence of lower eigenvectors). Assume that Ai(]HI) < infcress(m) 
for some L > 1. Then XlO^) is the L-th eigenvalue ofM and for N large 
enough, XLiH^) is the L-th eigenvalue of H^. Furthermore, if (^7v^)Ar>2 
is a sequence of associated eigenvectors, then, up to a subsequence. 



lim C/7v*S^^ = 



(21) 



strongly inF^, where ^^^^ is an eigenvector o/H associated with the eigen- 
value Xl{^. The latter convergence is strong in the norm induced by the 
quadratic form of M on if (A3s) holds true. 

The proof of Theorem [2] is provided in Section lY.ll Let us now make some 
comments on this result. 

The weak limit (|19p shows how the Bogoliubov Hamiltonian EI arises 
from the particle-conserving Hamiltonian H^. This convergence essentially 
implies the upper bounds on the eigenvalues of Hn and this is the easy part 
of our proof. In previous results (for instance in [55]) the upper bounds were 
more involved because the argument was carried out in the original Fock 
space J- and the trial state had to be projected on S^^ . On the contrary 
we work in the excited Fock space J-+ and only need to project the state 
on with M < iV (in practice 1 <C M < iV) before applying U^, 

which is much easier. Note that we actually do not need the condensation 
assumption (A3) for showing (i) and the upper bound on the eigenvalues. 

The difficult part of the proof is the lower bound on the eigenvalues, which 
requires (A3) as well as a localization method in the Fock space in the 
spirit of a previous result of Lieb and Solovej [4H Thm. A.l]. The idea is 
to estimate the error made by replacing a vector $ G by its truncation 
on J-^^ , in a lower bound. This method is the object of Section [6] where 
we prove an IMS-type localization formula. 

Even if our simplified statement does not reflect this, we are able to prove 
explicit error estimates. For instance, for the convergence of \l{His[) in (ii), 
we will prove that 

- CLisR^N) + N-^l^) <\l{Hn) - Nen- <ClN-^I^ (22) 

for N large enough, where Cl and Rl only depend on (Al(EI) — Ai(E[)) and 
{Xl{Hn) — Xi{Hn)), and where we recall that er^ is given in Assumption 
(A3). Similar estimates can be obtained for the eigenfunctions. Under the 
strong condensation assumption (A3s), our proof shows that one can take 
eji{N) = 0{N~^), leading to an overall error of the order 0(A^-V3). 

In [521l24j. Seiringer and Grech have shown a result similar to Theorem [21 
More precisely, in [52] Seiringer treated the case of O a cube in any dimen- 
sion, T = —A with periodic boundary conditions, and w a bounded and 
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positive periodic function such that w > 0. In ^24j, the same method was 
used to treat the case ofr2 = IR'^,T = — A + T4xt with Vext — )• oo at infinity, 
and w a bounded positive function such that w >0. In these two cases, the 
properties (Al), (A2) and (A3s) are all satisfied and Theorem [2] applies. The 
proof of strong condensation (A3s) is simple and relies on w being positive, 
see, e.g., [52l Lemma 1]. 

The speed of convergence ()22p that we can establish in our abstract setting 
is slightly worse than the 0{N-'^/'^) obtained in ^[M]. The method of 
proof is rather different, however. In [521124] . the authors relate to an A^- 
dependent Hamiltonian H^°^ which is quadratic in the effective annihilation 
operators bj = ajaQ/\/N — 1. These operators only satisfy the commutation 
relations in the limit N ^ oo. The effective Bogoliubov Hamiltonian H^°^ 
can be diagonalized by a unitary transform, up to an error. The unitary is 
constructed by inserting the effective 6j's in the formula of the Bogoliubov 
rotation which diagonalizes the true Hamiltonian H. 

In the present paper, by applying the unitary Un, we settle the problem in 
the excited Fock space in which the true Bogoliubov Hamiltonian H acts. 
The main advantage of our approach is that and H are now independent 
of A^, which makes the comparison with UnHnU'^ clearer, in our opinion. 
The effective operators 6j's were also used in previous works |4 H 1421 39j on 
the one- and two-component Bose gases, for which our approach could be 
useful as well. 

Let us remark that the convergence (|20p in the Fock space J-+ can be 



rewritten in the original A^-body space as 



N 



(1) 
N 







N 

E 



®(Af-i) 




where = </>q^'* © (p^^^ © • • • G In particular, when t;; / 0, we see that 
the many-body ground state ^'j^^ is never close to the Hartree state in 
the norm of This is because the vacuum in is never an eigenvector 
of the Bogoliubov Hamiltonian H, except when w = 0. A similar property 



holds for the lower excited states ^ 

,iL) 



(L) 
N 



Although the wavefunction ^ \t is in general not close to the Hartree 



state Uq in the norm of Sj , its density matrices are close to that of Uq . 
Indeed, if (Al)-(A2)-(A3s) hold, then the convergence of eigenvectors in (iv) 
implies the following convergence of the one-body density matrices 



lim Tr 

A^oo 



Qi^(L)Q 

^ N 



7$(-£.) 







and 



lim Tr 

A-i-co 



P7^(^)^-(A^-A^f^)l^^o)(^^o| 







where P = |tio)(no| = 1 — Q and = Tr7^{L). In particular, by simply 

controlling the cross term P^^^^Q by the Cauchy-Schwarz inequality, we 
immediately obtain the complete Bose-Einstein condensation 



Tr 



A^ 



-1 



t^o)(no| <0{N-y^). 
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More generally, we can show that if ^' is a wave function in S)^ satisfying 
(*7V,iJAf*Af) < E{N) + 0(1), then we have 



Tr 



T^^l-\uf){uf\ <0(iV-V2) (23) 



where A; G N is fixed, and rj^l^ is the operator on S)^ with kernel 



"ifNixi, ...,XN)'^N{yi, ■.■,yN)dxk+i...dxN dyk+i-.-dyN- 



Note that, by looking at the first order of the density matrices in (|23p . the 
excited states cannot be distinguished from the ground state. A slightly 
weaker version of (p3]) . namely that rj^^(xi, x^; xi, x^) converges to 
|no(xi)p...|no(xfc)p weakly in S)^, was proved recently by Kiessling for the 
ground state of bosonic atoms [30j. In fact, the convergence of reduced 
density matrices in ([23|) is well understood in the time-dependent setting [571 
[2m in[2l H9lll4j . and in this case there has been recent interest in corrections 
to the Hartree equation [26l [25l |6] , where our method might also apply. 

Remark 6. In the convergence in (in), the sequence of ground states {^tv^Iat 
can be replaced by any sequence of approximate ground states. More pre- 
cisely, if for every N we take a wave function "^'^ G such that 



lim UHn)^'-Xi{Hn)] =0, 

N--¥oo \ " / 

then we still have Un^'j^ — )• ^^^^ in the same sense as in Theorem\^ 

Remark 7. // we perturbate the factor 1/{N — 1) in front of the interaction 
term in the Hamiltonian by a term of 0{N''^) order, then our results 
in Theorem remain valid, with the modification that the eigenvalues are 
shifted by an extra term of order 0(1). More precisely, if kn = {N — 
V)~^ + kN~'^ + o{N'^'^) with K G M fixed, then the min-max values of the 
Hamiltonian 

N 

Hn,k := ^ Tj + kat ^ w{xi - Xj) 

i=l l<i<j<N 

satisfy 

XUHn,.) = Neu + K(/iH - en) + Al(BI) + o(l). 

2.5. Positive temperature. At a positive temperature /3~^ > 0, the free 
energy of the system is given by 



F^(iV) := inf { TT[HNr] - r'S{r)} = -^^ logTr^, 

Tr^,ivTr)=l 



where ^(r) := — Tr[riogr] stands for the von Neumann entropy. Our goal 
is to establish the convergence 

lim {Fp{N) - Nen) = -^^ logTr^+ 
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For this we need two new conditions. At positive temperature the stabihty 
of the system does not follow from the simple relative boundedness assump- 
tions (Al) on w. So we need the following 

(A4) (Stability). There exists (3o > such that F^o(^) ^ for all N 

and Ttjt^ [e-^o"^] < oo. 

Our second new assumption is a modified version of the zero-temperature 
condensation (A3), which we now only assume to hold for the Gibbs state 
at temperature I3~^ for simplicity. 

(A3') (Bose-Einstein condensation at positive temperature). For any /3~^ < 
, one has 

hm = 1 (24) 

where ^p^N is the one-body density matrix of the Gibbs state T^^^y ■= e~^^'^ 
/Tr [e~^^^], namely, in terms of kernels, 

ll3,N{x,y) := N j Ti3^N{x,X2,...,XN]y,X2,...,XN)dX2...dXN- 



Let us remark that if the strong condensation assumption (A3s) holds true 
for some Eq G (0, 1) and Trjr^ |^g-(i-eo)/3oHj ^ ^ some /3o, then we can 
prove (A3') and (A4) for the corresponding /Jq. We of course always assume 
that (Al) and (A2) hold true. Moreover, if /i > r/n and Ki = K2 > 0, then 

dr{h + C) + C > H > dr(/i -e)-Ce 

(see ()69p in Appendix[X|), and hence the condition Trjr^ |^g-(i-£o)/3oHj ^ ^ 
is equivalent to Tr |^e~^^~^°^'^"''] < 00. The latter holds true if we have 
Tr [e-(i-=")(i-°i)^o'^] < 00, because h>{l- ai)T-C, where qi G (0, 1) is 
given in Assumption (Al). 

Our main result is the following 

Theorem 3 (Positive temperature case). Assume that (A1)-(A2)-(A3')- 
(A4) hold true. Then for every I3~^ < , we have 



lim Trjr 

n— >oo 



N 



0. 



This implies the convergence of the corresponding Gibbs states and of the 
free energy: 



lim {F^{N) - Nen) = -^^ logTr^^ 

TV— >-oo 



Theorem [3] is proved using the same argument as that of the proof of The- 
orem [21 together with a well-known localization inequality for the entropy, 
see Section YT72\ 
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3. Applications 



It is easy to verify that the model considered by Seiringer in j52j and by 
Grech-Seiringer in [21] satisfy our assumptions (Al), (A2) and (A3s). So 
our approach apphes and we are able to recover their results. 

In the rest of this section, we consider two Coulomb systems for which we 
are able to apply our main results. 

3.1. Bosonic atoms. For a bosonic atom we mean a system including a 
classical nucleus at the origin in of charge Z > and "bosonic quantum 
electrons" of charge —1. The system is described by the Hamiltonian 



N 



-A,; 



z 

I 



+ 



E ■ 

l<i<j<N 



1 



acting on the symmetric space Sj'" = (g)^^^^L^(M^). For simplicity of writing, 
we only consider spinless electrons. We shall study the asymptotics of Hn^z 
when N ^ oo and (A^ — 1)/Z is fixed. 

By using the unitary Vi 
(^^/^^{Ixi, ...,£xn) with £ = N 



N 

E 

1=1 



-A,; 



Sj^" defined by {Vi^){xi, ...,xn) 
- 1, we can rescale the Hamiltonian to 

1 \ 1 1 



tlXi 



+ 



iV 



E 



1 ^ — ' \Xi 

l<i<j<N ' 



where t := (A^ — 1)/Z. This Hamiltonian has the same form as in the 
previous section, with Q, = M'^, T = —A — l/(t|x|) and w(x — y) = \x — y\~^. 
Note that Assumption (Al) holds due to Kato's inequality 



1 

— < 



eAx H on L 



4e 



(25) 



and the fact that Ht^N is self-adjoint on (^^^i7^(M^). Note also that the 
first eigenvalue of Ht^N, when it exists, is always non-degenerate. In the 

following we shall denote by ^[^li the corresponding unique positive ground 
state. 

It was already proved by Benguria and Lieb [7| that for every t > 0, the 
leading term of the ground state energy of Ht^N is given by Hartree's energy, 
that is 



mia{Ht,N) = Nen{t) + o{N) 



as — )• oo, where 



enit) : 
Here, again 



inf 

ll«ll=i 



I Vii(x) 



'dx 

t 



\x\ 2 



\u\ 



(26) 



.(27) 



Dif,9) :-- 



f{x)g{y) 
\x - y\ 



dx dy. 



Note that D{f,f) > because the Coulomb potential \x\ ^ has positive 
Fourier transform. 
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The fact that the next order in the expansion of the ground state energy of 
bosonic atoms is given by BogoHubov's theory was first conjectured in |44| . 
In the fohowing, by applying Theorem [2l we shall establish not only this 
conjecture but also many other properties of the system. 

By a convexity argument, it can be shown (see [34j ) that the Hartree 
minimization problem ()27p has a minimizer ut if and only if t < tc, for some 
critical number tc G (1,2) (it was numerically computed in [5j that tc ~ 
1.21). In the case of existence, the minimizer is unique, positive, radially- 
symmetric. Moreover it decays exponentially and it solves the mean-field 
equation 

ht Ut = 0, 

/li := -A - ^ + |nt|2 * ^ - /iH(t), 

with the Lagrange multiplier finit) < 0. Moreover, if t < tc, then fiH{t) < 
and there is a constant 7?H(i) > such that 

ht > m{t) > on i3+ := {ut}^. (28) 

The critical binding number tc in Hartree's theory also plays an important 
role for the original quantum problem. In fact, it was shown in [71|5ll[3] that 
for every there are two numbers b{N) < b'{N) satisfying that Ht^z always 
has a ground state if t < b{N) and Ht^z has no ground states if t > b'{N), 
and that 

lim b{N) = lim b\N) = tc- 

N^oo Af— s-oo 

In the following we shall always assume that t is fixed strictly below tc- 
In this case. Assumption (A2) holds true. In fact, due to Hardy's inequality 

^ < -A^ on L^{R^), (29) 



4 X 



the function Kt{x,y) := ut{x)\x — y\~^ut{y) belongs to L^((R^)^). Hence, 
Kt{x,y) is the integral kernel of a Hilbert-Schmidt operator, still denoted 
by Kt. Note that Kt > because \x — y\~^ is a positive kernel. Thus the 
spectral gap (|28|) implies the non-degeneracy of the Hessian, namely 



ht + Kt Kt 
Kt 



The condensation in Assumption (A3) is implicitly contained in the proof 
of the asymptotic formula (f26|) by Benguria and Lieb. In fact, the upper 
bound in (j26p can be seen easily by using the Hartree state uf^ . The 
lower bound is more involved and it follows from the Lieb-Oxford inequality 
, [37] which says that for every wave function ^' G S^^ , 



^' I ^ T W ) ^ Jl)(p*,p*) - 1.68 / p^ixf'^dx. (30) 



The following is a quantitative version of the strong condensation (A3s). 
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Lemma 4 (Strong condensation of bosonic atoms). If t <tc, then 

N 

Ht,N - Nen > (l - N~^/'') Y.^ht)i - CN'/\ (31) 

i=l 

Remark 8. In particular, from [31\) it follows that if (-fft,Ar)* < Nea + R, 
then (uo,7*Mo) > - CN^/^ which is (A3). When t < t^, by Theorem\M 
we can improve the estimate to (uo,7>iitio) > N — C. The latter was shown 
by Bach, lewis, Lieb and Siedentop in [A\ using a different method. 

Proof. We start by estimating the terms on the right side of the Lieb-Oxford 
inequahty (j30p . First, from the positivity D{f,f) > 0, we have 

\d{p^,P^) > D{p^,N\ut\^)-^D{N\ut\\N\ut\^) 

= ND{p^,\ut\^)+N\en{t)- finit)). (32) 

On the other hand, using the Hoffmann-Ostenhof inequahty [29| and Sobolev's 
inequahty [36\ Theorem 8.3] we can estimate 

Therefore, by Holder inequahty, we find that 

I pT<[j P%Y' (/ P*)'^'<eTr[-A7d + C^. 

Thus using the Lieb-Oxford inequahty ([50]) . the estimate ([52]) and ht > 
-A/2 - C we get 

{Ht,N)^ > Nenit) + (l " §) T^[ha^] " 2Ce - Cs-'N^^^ 

for all e >0. Replacing e by we obtain ([311). □ 

All this shows that if t < tc, then Assumptions (Al)-(A2)-(A3s) hold 
true, and we may apply Theorem [2] to show that the lower spectrum of 
converges to the lower spectrum of the Bogoliubov Hamiltonian 

a*{x){{ht + Kt)a) {x)dx 

+ ^ j^j^^^'^^'^^i^*^^^"*'^^^ ^ a(x)a(y)^ dxdy, 

which acts on the Fock space F+ = ®'^=Q®^yja^+- Beside some basic 
properties of already given in Theorem[Tl we have the following additional 
information. 

Proposition 5 (Bogoliubov Hamiltonian of bosonic atoms). For every t G 
(0, tc) one has 

aess(IHO = [infcj(Mt)-//H(t),oo). 

Moreover, ift < 1, thenMt has infinitely many eigenvalues below its essential 
spectrum. On the other hand, if t > 1, then Mt only has finitely many 
eigenvalues below its essential spectrum. 
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Proof. First, since ht = — A+Vt{x) — ^u(t) , where Vt{x) := {\uo\'^ *\.\^^){x) — 
t~^\x\~^ is relatively compact with respect to —A, we obtain (Tess(^i) = 
[— //h(O)Oo)- Thus cress(EIi) = [infcr(]HIt) — /XH(t),oo) by (ii) in Theorem[TJ 
The other statements follow from (iii) in Theorem [1] and the fact that ht has 
infinitely many eigenvalues below its essential spectrum if and only if t < 1. 
In fact, if t < 1, then Vt{x) = mt* \.\~^ where the measure 

has negative mass J m < 0, and we can follow [4'6\ Lemma II. 1]. 

On the other hand, if t > 1, then by applying Newton's Theorem for the 
radially symmetric function |no(x)p, we can write 

vd^) > I I";'","; - ^ = / i...(<,)i^ (if - i) d,. 

JiRS max{|x|, |y|} \x\ J\y\>\x\ \\y\ \x\J 

Because ut{x) decays exponentially, we obtain that also decays ex- 

ponentially and hence belongs to L^/^(M^). Therefore, by the CLR bound 
(see e.g. |38l Theorem 4.1]), we conclude that ht has finitely many eigenval- 
ues below — /iH(i)- n 



By the celebrated HVZ Theorem [47], one has 

(^ess{Ht,N) = p4,Ar,Oo) 

where St, at is the ground state energy of the (A^ — l)-body Hamiltonian 

Af-l 



1 \ 1 >p 1 



Xi 



=1 ^ ' l<i<j<N-l 

An asymptotic formula for T,t^N can be obtained by the same method as 
that of Theorem [2] (see Remark [7]), namely 

St,7v = Nenit) - mit) + inf cT(Mt) + o(l) = Neuit) + mi a,U^t) + o{l). 

It thus turns out that the essential spectrum of Ht^N — ^euit) converges to 
the essential spectrum of Mt- 

Thus, by combining the HVZ Theorem and Theorem [21 we obtain the 
convergence of the whole spectrum of Ht^N as follows. 

Theorem 6 (Spectrum of bosonic atoms). If t £ (0,tc) is fixed, then we 
have the following statements in the limit N ^ oo. 
(i) (Essential spectrum) . The essential spectra of Ht^N cLnd are 
aess{Ht,N) = [St,Ar,oo), o-ess(IHIt) = [inf (T(IHIt) - fiu{t),oo) 
and we have the convergence 

lim {^t,N - Nen{t)) = miaiUt) - m(.t). 

(a) (Ground state energy and ground state). We have 

lim {ini a{HtN) - Neu{t)) = inf (T(IHt) < 0. 

TV— >oo ' 

For N large enough, Ht^N has a ground state ^'j^^f (which is then unique 
and can he assumed positive ) and 

lim C/^^S^i = ci>W 
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strongly in the norm induced by the kinetic energy in Tj^, where <I>j"^^ G J> is 

the unique ground state o/H such that > 0, Q, being the vacuum 

state in Fj^ . 

(Hi) (Lower eigenvalues and eigenstates) . Assume that Mf has L eigenvalues 
below its essential spectrum, for some L > 1. Then for N large enough, Ht^N 
also has L eigenvalues below its essential spectrum, and the L-th eigenvalue 
of {Ht^N — Ne-B_{t)) converges to the L-th eigenvalue of M.t as N ^ oo. 
Moreover, if ^'j^ is a corresponding eigenvector, then up to a subsequence. 



Uiy^i^}, converges to an eigenvector of the L-th eigenvalue ofMt as N 



oo. 



3.2. Trapped Coulomb gases. We consider a two- or three- dimensional 
system of N bosons trapped by an external potential and interacting via the 
Coulomb potential. The system is described by the Hamiltonian 



N 

wixi 



acting on the symmetric space Sj^ = L^(M'^) for d = 2 or d = 3, where 

/C = — A or /C = Vl — — 1 is the kinetic energy operator. 

Our study in this section was motivated by a recent work of Sandier and 
Serfaty |5H I50j . They proved in [50] that for the classical 2D Coulomb 
systems, the ground state energy behaves like 

Nen -jlogN + c + o{l) 

where en is the Hartree ground state energy when /C = and c is the ground 
state energy of an infinite system of point charges in the plane, interacting 
through the Coulomb potential and with a neutralizing jellium background. 
The log term is due to some local scaling invariance (think of replacing 
y by a constant in a small box and shrinking the system about the center 
of this box). We shall see that, when the kinetic energy is introduced, since 
the local scaling invariance is broken and shrinking the particles has a too 
large cost, the logarithmic term disappears and the next term is of order 
0(1), given by Bogoliubov's theory. In order to make this vague statement 
rigorous, we need some natural assumptions on T = fC -\-V and on w. 



Assumptions. Ifd = 2, then w{x) = — log /C is either —A or \/—A -\- 1- 

loc( 



1, and V G L? (M^,M) satisfying that 



liminf - — V4pj > 0. 

I^I^oo [log |a;|J^ 

If d = 2>, then w{x) = \x\~^ , K, is either —A or \/—A + 1 — 1, and V £ 
Lf^^iM.^ ,W) satisfying that V{x) oo as |x| — )• oo. 

Let us show that under these conditions. Assumptions (Al)-(A2)-(A3s) 
in Theorem [2] hold true. First, note that }C -\-V is bounded from below on 
L2(R'^) (see [361 Sec. 11.3]). When d = 3, K = - 1, then the 

relative bound ([5]) in (Al) follows from the relativistic Hardy inequality \38\ 
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Lemma 8.2] 

<V^< VT^^ on L^CR^), 

7r|x| 

while if /C = — A, then we even have the stronger bound 

\w{x-y)\^ <Co[}C^ + }Cy + V{x) + V{y)] + C on L^iR'^f) (33) 

due to Hardy's inequahty When d = 2, then ([55]) also holds true, due 
to the estimates IC + V > C-i[log(l + \x\)]^ - C and 



\w{x-y)\'^ < 2 



^ +[log(l + |x|)]2 + [log(l + |y|)]2 



\x - y\ 



+ C 



for some C > 0. Thus (Al) holds true. On the other hand, Assumption 
(A2) follows from the following 

Proposition 7 (Hartree theory). Under the previous assumptions, the vari- 
ational problem 

en:= inf \ {u, {)C + V)u) + ^ [ [ \u{x)\Mx - y)\u{y)\^ dx dy] 
IMI=i 

has a unique minimizer uq which satisfies that uq{x) > for a.e. x G M*^ 
and solves the mean-field equation 

huQ = 0, 

h := IC + V + \uq\'^ * w — nu, 

for some Lagrange multiplier /hh G R- The operator h has only discrete 
spectrum Xi{h) < X2ih) < Xsih) < ... with limj_j,oo Aj (/i) = oo. Moreover, 
the operator K with kernel K{x,y) = uo{x)w{x — y)uo{y) is Hilbert- Schmidt 
on L^(M'^) and it is positive on S^^. Finally, uq is non- degenerate in the 
sense of ([7]). 

Before proving Proposition \7\ let us mention that in 2D, the Coulomb 
potential w{x) = — log|x| does not have positive Fourier transform. More 
precisely, w = pv \ ■ \~^, the principal value of | • |~^. Although vu is not 
a positive type kernel, we still have the following restricted positivity (see 

m). 



Proposition 8 (Coulomb log kernel). For any function f G L^{R'^) D 
L^^'^{R'^) for some e > with 

/ log(2 + |x|)|/(x)| dx < oo and / / = 0, 



we have 



< D{f,f) ■=- jj f{x)log\x-y\f{y)dxdy < oo. 
We can now provide the 
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Proof of Proposition^ By using the relative bound ([5]), Sobolev's embed- 
ding [361 Theorem 8.4] and the fact that V{x) — )• oo as |x| — )• oo, it is 
straightforward to show that there is a minimizer uq for en in (j34p such that 
uo G L>(/Ci/2) and 

oo. 

Since {u,ICu) > {\u\,IC\u\) [361 Theorems 7.8 and 7.13], we may assume that 
Uq > 0. Moreover, the Hartree equation huo = imphes that uo{x) > for 
a.e. X G and no is then the unique ground state of h. Moreover, because 
V{x) + (luoP * w){x) — 7- oo as |x| — >• oo, the operator h (on has only 
discrete spectrum Xi{h) < \2{h) < ... with \i{h) > and limj_j.oo = oo 
(see [m Theorem XIII. 16]). 

For any normalized function u £ S), we have J(|tip — |tioP) = and hence 
D{\u\'^ — |noP, |np — |noP) > by Proposition [SI Using this and a convexity 
argument we can show that uq is the unique minimizer of en in (j34p . 

The operator K with kernel K{x,y) = uo{x)w{x — y)uo{y) is Hilbert- 
Schmidt in L'^{{R'^f). When d = 2, or d = 3 and /C = -A, this fact 
holds true due to ([M]). When d = 3 and /C = Vl - A - 1, the Hilbert- 
Schmidt property follows from the Hardy-Littlewood-Sobolev inequality [361 
Theorem 4.3] and the Sobolev's embedding H'^/'^{R^) C L3(m3). 

Finally, the operator K is positive on S)^ because {v, Kv) = D{vuq, vuq) > 
for every v G The latter inequality follows from Proposition [8] and the 
fact that / (vuo) = 0. We then deduce from h > \2{h) — \i{h) := r/H > 
that Uq is non-degenerate in the sense of ([7]). □ 

We have shown that (Al) and (A2) hold true. Therefore, we may consider 
the Bogoliubov Hamiltonian 

n=j a*(x)(/i + K)a(x)dx + ^ jj K{x,y)(a*{x)a*{y) + a{x)a{y)^dxdy, 

which acts on the Fock space = '2)"ym-^+- From Theorem [1] and 

the spectral property of h, the following can easily be proved. 

Proposition 9 (Bogoliubov Hamiltonian of trapped Coulomb gases). Un- 
der the above assumptions on /C, V and w, the Bogoliubov Hamiltonian M. is 
bounded from below. Its spectrum is purely discrete, consisting of a sequence 
of eigenvalues Ai(E[) < A2(]HI) < A3(EI) < • • • with \\m.j^oo Aj(EI) = oo. 

Now we consider Assumption (A3s). In three dimensions, the condensa- 
tion can be obtained by following the proof of LemmaO In two dimensions, 
we have the following result. 

Lemma 10 (Strong condensation of 2D Coulomb gases). Under the above 
assumptions on K., V and w, we have 

N 
i=l 

This condensation can be proved by using the same argument as in the 
proof of Lemma H| and the following Lieb-Oxford type inequality, whose 
proof can be found in Appendix [Q 
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Proposition 11 (Logarithmic Lieb-Oxford inequality). For any wave func- 
tion ^ G L^(]R^) such that l^'p is symmetric and £ n for 
some < e < 1 and f^2 log(2 + \x\)p^{x)dx < oo, we have 




Xj\j > ^D{p^,p^) - ^ J^^p^logp^ 
n / r \ r / n 



e 



C p^-- p^ . (35) 



/k2 



The estimate (|35p is probably not optimal but it is sufficient for our 
purpose. In particular we do not know if the error term involving e can 
be removed. By applying Theorem [2] and Theorem [3] we get the following 
results. 

Theorem 12 (Spectrum of trapped Coulomb gases). Assume that IC, V 
and w satisfy the above assumptions. 

{i) (Eigenvalues). The Hamiltonian H]\j has only discrete spectrum Xi{Hf^) < 
X2{Hn) < X3{H]\f) < ... with limL^oo ^l{Hn) = oo. For every L = 1,2,... 
we have the following convergence 

lim {Xl{Hn) - Neu) = 

where Xl{M) is the L-th eigenvalue of the Bogoliuhov Hamiltonian H. 

{ii) (Eigenvectors). If '^^^^ is an eigenvector of the L-th eigenvalue of H^, 

then up to a subsequence, Un'^^^^ converges (strongly in the norm induced 
by KL + V ) to an eigenvector of the L-th eigenvalue o/Ht as N ^ oo, where 
Un is defined in [T^) . 

(Hi) (Positive temperature) . If we assume furthermore that Tr e~^'^'^'^^^^ < 
oo for some /3o > 0, then for every < < I3q^ we obtain the convergence 



lim Trr, 

N^oo ^ 



= 

In particular, we have the convergence of the free energy 



lim ( - log Tr^,^ e"''^^ - N en) = -T^ log Tr^, e'^^. (36) 



We have a couple of remarks on Theorem [T5J 

First, we note that the condition Tr e"'^"^'^"''^-' < oo is satisfied if V grows 
fast enough at infinity. For example, if/C = — A,(i = 2ord = 3, and 

r . n y{x) d 
hm mf j — - > — , 

|x|-i.oo log |x| po 

then one has the Golden-Thompson-Symanzik inequality l63| [61] (see 
also [18] for an elementary proof) 



Tre-*('C+^)<(47r/3o)-'^/' / e-^^^'^^U 



Ix < oo. 



Moreover, if /C = a/— A + 1 — 1, d = 2 or c? = 3, and 

limmf > 0, 
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then Tie-'^oiic+v) ^ ^ f^^, > g, due to [181 Theorem 1] and the 

operator inequahty \/—A + \x\ > y^— A + |xP in L^(M'^). The latter is a 
consequence of the operator monotonicity of the square root and the fact 
that V^\x\ + \x\y/^ > in L^iW^), see JW, Theorem 1]. 

Second, in [50l Theorem 1] the authors provided upper and lower bounds 
on the classical free energy at nonzero temperatures, which coincide when 
I3~^ — )• 0. In the quantum case we are able to identify precisely the hmit (j36p 
even when /3~^ > 0, which is given by the Bogohubov Hamiltonian. 



4. Operators on Fock spaces 

In this prehminary section, we introduce some useful operators on Fock 
spaces and we consider the unitary defined in (jlSp in detail. 

For any vector / G i^, we may define the annihilation operator a{f) and 
the creation operator a*{f) on the Fock space T = 0^=o^^ follow- 
ing actions 

a{f) j Yl /-(D^-^Z-W I = ^ E (/./<x(i)>/<x(2)®-®/<x(iV),(37) 
\cre6jv / o-e&N 



«*(/^) Yl Ml) ® - ® faiN-1) I = -^-(1) ® - ® (38) 

\crGejv_i 

for all /, /i, /at in -5, and all A'" = 0, 1, 2, .... These operators satisfy the 
canonical commutation relations 

[a{f),aig)]=0, [a*if),a*ig)] = 0, [aif),a*ig)] = {f,g)n. (39) 

Note that when / € Sj-^, then a(/) and «*(/) leave J-+ invariant, and 
hence we use the same notations for annihilation and creation operators on 
J^f . The operator- valued distributions o(x) and a*{x) we have used in ([12]) 
can be defined so that for all / € 



a{f) = / f{x)a{x)dx and a*{f) = / f{x)a*{x)dx. 
Jn Jn 

To simplify the notation, let us denote a„ = a{un) and a* = a*{un), where 
{un}^=0 is an orthonormal basis for L^(il) such that uq is the Hartree min- 
imizer and Un G D{h) for every n = 1,2,.... Then the Bogoliubov Hamil- 
tonian defined in (1121) can be rewritten as 



{Um,{h + Ki)Un)L2(^Q)a^an + ^ (Urn's) Un,K2)i^2(^Q2) al^On 

m,n>l 

+ ^{K2,Um'SUn)L^(n-^)aman- (40) 

The sums here are not convergent in the operator sense. They are well de- 
fined as quadratic forms on the domain given in (jl4p . Since the so-obtained 
operator is bounded from below (by Theorem [1]), it can then be properly 
defined as a self-adjoint operator by the Friedrichs extension. 
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It is also useful to lift operators on S)'^ to the Fock space T. The following 
identities are well-known; their proofs are elementary and can be found, e.g., 
in [8] and [56l Lemmas 7.8 and 7.12]. 

Lemma 13 (Second quantizations of one- and two-body operators). Let 
A be a symmetric operator on S) such that Un G D{A) for all n > 0, and 
let uj be a symmetric operator on Sj such that Um ® Un ^ D[uj) and 
{um CS) Un, uUp® Uq) = {un fX" Um, UJ Up Ug) for all m, n,p,q > 0. Then 

oo N 

dT{A) := Oe = {um,Aun)sj a*^an 

N=l j=l m,n>0 

and 

oo ^ 

^ ^ 2 ^^"^ n„, a; Up Ug)^y2 al^a^apaq 

N=21<i<j<N ■m,n,p,q>0 

as quadratic forms on the domain 

oo M N 

U 0(2)(Span{no,ni,...}) C-F. 

M=0 N=0 sym 

The same identities also hold for operators on the Fock space where 
we can use the orthonormal basis {un}^=i for In particular, the particle 
number operator M = dr(l) = YlJLoj^Sji ^ can be rewritten as M = 
Yl'^=o 'Ki'^n, and the particle number operator on F+ is 7V+ = Yl^=i 'Ki'^n- 

By using the second quantization, we can write Hj^ : Sj^ — )• S^^ as 

^ \ ^ TmnamCLn + ^, ^ _ ^, ^ Wmnp<?aman«p«g I • (41) 

1 m,n>0 m,n,p,q>0 

Here Tmn ■= {um,Tun) and 



Wmnpq ■■= // Um{x)Un{y)w{x - y)Up{x)Uq{y)dxdy . 

Since Hn and EI live in different Hilbert spaces, to compare them we 
need to use the unitary transformation Un : Sj^ — )• J-^^ defined before in 
(|18|) . The action of Un on annihilation and creation operators is given in 
the following lemma, whose proof is elementary and is left to the reader. 

Proposition 14 (Properties of Un). The operator Un defined in (jlSp can 
be equivalently written as 

N / N~j \ 

f/^W = 0Q^^(^^^^^j (42) 
for all G S)^ , and where Q = I — \uo){uq\. Similarly we have 
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'J 



■T-<N 

on J-T 



for all (f>j G j = 0, ...,N. These operators satisfy the following identities 

UNalaoUli = N-M+, 
UNa*{f)aoU*^ = a*{f)^N-M+, 
UNala{f)Ut, = VN-Af+a{f), 
UNa*if)aig)U*^ = a*{f)aig), 

for all f,g £ S)+. 

Note that the previous properties are purely algebraic. They do not 
depend on any special choice of the reference one-body function uq £ S). 
Roughly speaking, the unitary transformation Un{-)U^ leaves any or 
am invariant when m > 1, and it replaces each oq, and Og, by ^/N — N+. 
The later is essentially the number \/N if Af+ is small in comparison with 
A^, on the considered sequence of states. 

By using these identities and the commutation relations (j39p , we can com- 
pute easily Uj^AU^j where A is any operator on 9)^ which is a (particle- 
conserving) polynomial in the creation and annihilation operators. For ex- 
ample, 

UN{alalaman)U*N = UN{alajn)U'^UN{alan)U^ 
Using that a{f)^/N -J\I+ = y/N -M+ - la(/) for all / G S)+, we obtain 



A tedious but straightforward computation shows that 



4 



UNHMUtj-Nen = Y,^3 (44) 



where 



1 AA+(AA+ - 1) 
^0 - 2*^0000 _i 



•m 



N -M+-1 



N - 1 
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^2 = ^ {um, {T - fin) Un) a*^an 

m,n>l 

N -M+ 



+ ^ {Um, {\uo\'^ * W + Ki) 



N -I 

m,n>l 



+ 2 2^ {um®Un,K2)a^a^ 

m,n>l 



, Ura ® tin) , ^ ^ 

m,n>l 

m,n,p>l 



N _ 

m,?i,p>l 

^4 = 2{N — \) ^mnpqlrn'^n'^p'^q- 

m,n,p,q>l 

Here recall that en = Too + (l/2)l^oooo and /in = Too + Woooo- 

In the next section, we shall carefully estimate all the terms of the right 
side of (HUl to show that 



UnHnUn - Neu w H 

in the regime A/^ <C N. 

5. Bound on truncated Fock space 

The main result in this section is the following bound. 

Proposition 15 (Preliminary bound on truncated Fock space). Assume 
that (Al) and (A2) hold true. For any vector^ in the quadratic form domain 
ofM such that $ S J^^ for some 1 < M < N , we have 

fW 

{UN{HN-Nenm)^-{m)^ <C^—{M + C)^. 

Here for a self adjoint operator A and an element a in a Hilbert space, 
we write {A)a instead of {a,Aa) for short. Recall that the quadratic form 
domain of H is the same as that of dT{h + 1), on which the quadratic form 
UN{HN—Neu)U^ is well-defined. As an easy consequence of Proposition [151 
we can prove the weak convergence in the first statement of Theorem [2j 

Corollary 16 (Weak convergence towards M). Assume that (Al) and (A2) 
hold true. Then we have for all fixed in the quadratic form domain of 
the Bogoliubov Hamiltonian H, 

lim {<^',UN{HN-Neu)Ut,^)^ =($',]HI$>^ (45) 
where by convention is extended to outside of J-^^ . 
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Proof of Corollary \ 161 It suffices to show the statement for <!>' = <I>. Note 
that the quadratic form domain of EI is the same as that of dT{h+l) because 
of (jl6p . and dT{h + 1) preserves all the subspaces iolp's. Therefore, if we 
denote by the projection of $ onto Tf'^^ , then 

lim = (]H)$ and lim lim (M + C)$^_$^, = 0. 

M— >-oo Af->oo Af— >oo 

If we denote Hn '■= Un{Hn — ^eH)U^, then from Proposition [T5] we 
have 

lim {Hjy — ]HI)$ = and lim lim {Hj\[)^ = 0. 

M — s>oo A/— >-oo Af— s>oo 

The latter convergence still holds true with Hn replaced by a non-negative 
operator H'j^ := ffjv + Co(]HI + Cq), where Co > is chosen large enough. By 
using the Cauchy-Schwarz inequality for the operator H'j^, we deduce that 

lim lim (i?Ar)$ - (-ffAf)$M = 0- 
M -^oo N—^oo 

We then can conclude that {Hjy),^ — )• (1HI)$ as — )• oo. □ 

The rest of the section is devoted to the proof of Proposition I15[ We shall 
need the following technical result. 

Lemma 17. // (A1)-(A2) hold true, then we have the operator inequalities 
on Fj^ : 

dr{QTQ) < — M + + C, 
1 — ai 

dr{Q{\uo\^ * \w\)Q) < -^^M + CM+ + C, 

1 — ai 

where 1 > ai > and a2 > are given in the relative bound ^ in Assump- 
tion (AI). 

Proof Using ([5]) we get |noP * w > —ai{T + C) and |tioP * \w\ < a2{T + C). 
Consequently, T < (1 — ai)^^h + C and |noP * \w\ < 02(1 — ai)~^h + C. 
The desired estimates follows from the lower bound IHI > dr(/i) — CA/V — C 
(see Remark nop . □ 

Now we give the 

Proof of Proposition \15[ Let $ be a normalized vector in the quadratic form 
domain of H such that ^ E J^^*^ for some 1 < M < N. Starting from 
the identity (jM]), we shall compare (^2)<i> with (11)$, and show that the 
remaining part is negligible. 

Step 1. Main part of the Hamiltonian. 
Lemma 18 (Bound on — H). We have 

(^2)$ - (M)$ < + c{N+ + 1)$ 

iV — i y i — ai 
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Proof. From (j40p and (j44p we have 

|2 



29 



dr (Q(|nor *u')Q + i^: 



AA, - 1 



iV- 1 

+5}? ^ (Um 1^ Wn,i^2) (amfln^)* 



(46) 



m,n>l 



where X := ^(iV - AA+)(7V - AA+ - l)/(iV - 1) - 1. 

Since dr (Qd^oP * w)Q + Kij commutes with J\f+, we have 

< dr {Qi\uo\^ * w)Q + K,) ^^^fi < ^^dr (Qdnol' * w)Q + K^) 

on Tt'^^ . By using Lemma [T71 and the boundedness of Ki, we can bound 



the first term of the right side of ()46p as 



dr {Q{\uo\^*nj)Q + K,)^^^±—^ 



< 



M 



N -I V 1 - ai 



02 



C 



(47) 



The second term of the right side of (|46p can be estimated using the 
Cauchy-Schwarz inequahty 



m,n>l 

1/2 



1/2 



< 



< 



< 



^ (g) Un,K2)Y 



2\l/2 



, m,n>l 



, m,n>l 



CM 



(48) 



Here we have used the fact that the operator K2 is Hilbert-Schmidt and the 
inequahty (X - 1)^ < 4(A/'+ + lf/{N - if and the estimate N+ < M on 



J"!*^ From dM]), (gZl) and dm), the bound in Lemma [H] follows. 



□ 



Step 2. Unimportant parts of the Hamiltonian. Now we estimate the other 
terms of the right side of (j44p . First at all, by using Ml_ < MM-\- on T^^^ 
we get 

(AA+(AA+-1))„ 



0000" 



< Af-i.. 



N -1 

The terms (^4)$ and (^3)$ are treated in the next lemmas. 

Lemma 19 (Bound on Ai). We have 

\{A,)^\<C^{M+)^. 



(49) 
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Proof. By using Hartree's equation (T + \uq\'^ * w — /Uh) uq = we obtain 

TmO = {Um,TuQ) = - (Um, (|uoP * w)uo) = - WmmO 

for all m > 1. Therefore, 

m>l 

By using the CauchySchwarz inequality we get 

1/2/ \ 1/2 



\m>l j \m>l 



Here in the last estimate we have used the bound 

2 



i^mooop — 



m>l m>l 

< 




n Jn 

2 

\uov-^)\^\uo{y)\^w{x - yf( 

n Jn 



< OO 



and the inequality Ml{N - M+) < M NM+ on T^''' . D 
Lemma 20 (Bound on A4). We have 

I (^4)$ I < f (M)$ + C(AA+)$ + C 

' ' iV — 1 y 1 — ai 

Proof. From Assumption (Al) we have 

Q^Q {a2(T^ + Ty + C)-w{x-y))Q^Q>0 

on i^^, where Q = 1 — \uq){uo\. By taking the second quantization (see 
Lemma [TH|) of this two-body operator, we can bound A4 from above by 

2(N- 1) ^ {um®Un,a2{T ® 1 + l(g)T + C)up0Uq) a'^a':^apaq 

m,n,p,q>l 

-dTiQTQ)iAf+ - 1) + l77T^-T^f+i^f+ - 1). 



N -1 ^ ' ' 2{N -I) 

By using Lemma [T7] and the fact that dT{QTQ) commutes with A/V, we 
obtain the operator inequality 

M ( "2 , ^AT , T-<J\/ 



N — 1 \1 — ai J 

Employing this argument with w replaced by —w, we obtain the same upper 
bound with A4 replaced by —^4, and the bound in Lemma [20] follows. □ 

Lemma 21 (Bound on A3). We have 
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Proof. Let us write w = W-^- — W- where w+ = max{w, 0} and W- = 
max{— WjO}. By taking the second quantization (see Lemma [T3]) of the 
non-negative two-body operator 



QCd{Q- eP)w+Q ®{Q- eP) + {Q - eP) ® Qw+{Q -eP)®Q 
+ {Q(^{Q + eP)w-Q (Q + eP) + (Q + eP) Qw-{Q + eP) ^ Q 



where Q = l — P = l — \uo){uo\ and e > 0, we obtain the following Cauchy- 
Schwarz inequality 

{um ^ Un,wup (g) uo)a^alapao 

m,n,p>l 

+ ^ (uo ® Up, wun ® Um)aoa*anam 

m,n,p>l 



{Um (K) Un, \W\Up Ugja^anapaq 

m,,n,p,q>l 



+e E {um O no, \w\un O no)amaoa„ao. 

m,n>l 

After performing the unitary transformation [/jv(-)t^Ar we find that 

^3 < {Um®Un,\w\Up®Uq)a*^a*^apaq 
m,n,p,q>l 

^ dT{Q{\uo\^*\w\)Q){N (50) 



N - 1 



for every e > 0. 

On the other hand, by the same proof of Lemma [20] we have 

{um®Un,\w\up®Uq){ama*^apaq)^ (51) 

m,n,p,g>l 

Moreover, by using Lemma [T71 we get 

dr(Q(|noP * \w\)Q){N - M+) <{N- l)dr(Q(|noP * \w\)Q) 
and , we get 

^dr(Q(|no|2*|H)Q)^^^^)^ < (dr(Q(|noP*|H)Q)>$ (53) 

VI - ai 

From ([50]) . ([5T]) and ([53]) . we can deduce that 



By repeating the above proof with w replaced by —w, we obtain the same 
upper bound on —(^3)$ and then finish the proof of Lemma [211 □ 
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Step 3. Conclusion. Using Lemmas [T8l [T9l [20l [2T] and ()i9]) . we can conclude 
from (Hil) that 



where ai and 02 are the constants appearing in the relative bound in (Al). 
Since Af+ < C{M+C) due to p^ . the estimate in PropositionfTSlfollows. □ 

6. Localization in Fock space 
From Proposition [TCI we have 

lim {UnHnU^ - iVen - H)$ = 

N^oo 

for every $ G J-^^^ with M fixed. In the next step, we want to localize a 
state in J^p into the smaller truncated Fock space J"^'" with M < N, 
without changing the energy too much. The following result is an adaption 
of the localization method used by Lieb and Solovej in [41J Theorem A.l]. 

Proposition 22 (Localization of band operators in J-'+). Let A be a non- 
negative operator on T such that PjD(A) C D{A) and PiAPj = when 
\i — j\ is larger than some constant a, where Pj is the projection onto i^-^. 
Let < f, g < 1 be smooth, real functions such that f'^ + g'^ = 1, f{x) = 1 
for \x\ < 1/2 and f[x) = for \x\ > 1 and let and gM be the localization 
operators on Tj^ defined by 

00 00 
fM = f{M+/M) = f{j/M)Pj and g.j = g{M+/M) = 9{j/M)Pj . 

j=0 j=0 

(i) We have 

^ ^0 < ^ - fuAfM - guAgM < -TTT-^o (54) 



where Aq = J^^j^^^j diagonal part of A andCj = 4(||/'||oo + ||5' 



00 1 



{a) Let Y be a finite- dimensional subspace of D{A) such that 
(dimy)~^ for every normalized vector <I> G y, then dim(/A/y) = dimY. 

The reader should really think of the discrete Laplacian (a = 1) in which 
case the inequality is nothing but a discrete version of the IMS formula 
[16[ Theorem 3.2]. The statement on the dimension of the localized space 
/a/ y will be useful to control eigenvalues via the min-max formula. A proof 
of Proposition [22] can be found in Appendix [Bj In the following we have two 
applications of this result in our particular situation. 

Our first application is a localization for H. 

Lemma 23 (Localization for H). For every 1 < M < N we have 
-^{M + C)<M- /mM/m - gM^OM <^{M + C). 

Proof We can apply Proposition [22] with ^ = EI - Ai(IHI) and a = 2. Note 
that the diagonal part of H is nothing but dT{h + Ki) which can be bounded 
from above by C(IH + C) , due to (HSj) . □ 
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Now we turn to a localization for '■= Un{Hn — Ne-ii)U'^. 

Lemma 24 (Localization for H]\f). For every 1 < M < N we have 
C 

' CdT{h): <N + C — Xi{Hn) ) < Hn — fuHNfu — guHNgu < 



M2 V ^ ^1-^^ 



<^(cdr(/i)|_^l^ + c7-Ai(F^; 



Proof. We apply Proposition [22] with A = H]\f — Xi{H]\f) and a = 2. From 
the inequality H]\r < C{M. + C) due to Proposition [151 we see that the 
diagonal part of Hj\f can be bounded from above by CdT(h).^<N + C. □ 

Remark 9. The error term dT{h)^j~<N in Lemma^2^ can be further replaced 
by C(Il + C) due to [W\). or replaced by C{Hn + CN) due to the stability 

N N 

HN>{l-ai)^Ti-CN>^^^hi-CN on i^^, 

^ l + a2^ 

1=1 1=1 



where ai and a2 are given in Assumption (Al). The bound dT{h) < C{H]^+ 
CN) on J-^^ also implies that Xi{H]\j) > —CN, although it is not optimal 
(we shall see that Xi{Hn) is of order 0{1)). 



7. Proof of Main Theorems 
We recall that we have introduced the notation Hj\f := U]\f{H]\f—Ne}i)U'^. 

7.1. Proof of Theorem [2]. The first statement of Theorem [2] was already 
proved in Corollary [16] above. We now turn to the proof of the other state- 
ments, using the localization method. 

Step 1. Convergence of min-max values. 

Upper bound. By applying Lemma[23|and Proposition I15[ we have, for every 
normalized vector <I> G D{Ml) and for every 1 < M < A^, 

H > /a/IH/m + 5m1H<7m-^(IH + C7) 



> fM 



fM 



+Xi{M)gl,-^{M + C). (55) 

Now fix an arbitrary L G N. For every 6 > small, we can find an 
L-dimensional subspace Y C D(IHI) C J-+ such that 

max < Xl(M) + 5. 

$gy,||$||=i 

By using the inequality gl^ < 2M+/M and 7V+ < C(IHI + C), which is due to 
(jl6p . we get ||g(jv/*^*|P 1^ Cl/M for every normalized vector ^ £Y. Here Cl 
is a constant depending only on Xl(E[) — Ai(EI). Therefore, if M > LCl, 
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then dim(/j\/y) = L by Proposition [22] (ii). Consequently, by the min-max 
principle, 

{fM^,HNfM^) 



max 

$ev||$||=i 



Thus from ()55p we obtain, for M large enough 

-1 



N 



Al(IH) > 1 



Cl 
M 



1 + Ci 



Cl 
M' 



By choosing M = N^^^ we get the upper bound 

Xl{Hn) - Nen = \l{Hn) < Al(H) + ClN'^I\ 



(56) 



Lower hound. By applying Lemma [23] and Proposition [T5l for every 1 < 
M < N we have the operator inequalities on 



<N 



C 



Hn > fuHNfu + guHNQM — + C — Xi{Hn, 



> I 



M 



1-C 




TM 



c 



+\i{HM)gli - ^(Cdr(/i)|^<iv + C-\i{Hm) 



(57) 



We first show that \i{Hn) > Ai(]HI) + o(l). For every 5 > small, we can 



find a normalized vector ^ D{Hn) C such that 
From Assumption (A3) we have 



N 



M 



M 



where er is as in (A3). We can choose R = 1. In particular, if we choose 



M > NeniN), then \ \gM'^ 



N\ 



as — )■ oo independently of 5 and the 



choice of ^N- Thus from ([FT]) and the simple bound dT{h) 



C{Hn + CN) (see Remark [9]) we can conclude that 

iV 
CN 



N\ 



^-c\r-]xi{m-c^^ 



M2 



+Xl{HN)\\gM^N\ 



By choosing M such that max{A^eR(A^), ViV} < M < A^, then taking 
(5 — )• 0, we obtain the lower bound Xi{Hn) > Ai(IHI) + o(l). 

By adapting the above argument, we can show that Xl^Hn) > Ai(EI) + 
o(l) for an arbitrary L G N. In fact, for every 6 > small, we can find an 
L-dimensional subspace Y C D{Hn) C such that 

max {HN)<i> < Xl(Hn) + 
$Gy,||*||=i 
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From Assumption (A3) and the upper bound Xl{Hj\[) — Xi{Hj\i) < Al(E[) 
Ai(EI) +o(l), we have 



ll5Af$|P < 



< 



M - M 

for every normalized vector ^ £ Y, where ept is as in (A3). We can choose 
R = Al(H) - Ai(]H) + 1. In particular, if M > NeR{N), then ||5m$|P ^ 
as A'^ — >■ oo independently of e and the choice of ^ in Y. Consequently, when 
N is large enough we have dim{ JmY) = L hy Proposition [22] (ii), and by 
the min-max principle, we get 



max 

$Gy,||*||=i 



mm' 

Thus from (l57|l and the simple bound dT(h),-i-< 
(see Remark [9|) we get 



> Al(BI). 

MHn) < C{Hn + CN) 



Xl{Hn) >{1-C 




Al(M) - C 




Cl 



M 



Cl 



N 

M2' 



(58) 



By choosing M such that vtiSiyi{N e r{N) , ^/ N] <^ M <^ N , we obtain 
Al(57v) > Al(IH) + o(1). 

Remark on the convergence rate. The error obtained in the lower bound ()58p 
is not better than \J £r{N) + A^~^/^. However, it can be improved by the 
following bootstrap argument. First, from ()57p with M = rN for some 
small fixed number r > 0, and the simple bound dr(/i)|T-<]v — Xi(Hn) < 

C{Hn + CN) (see Remark [9]), we obtain 

frNWrN<C{HN + C). (59) 

Next, by projecting the inequality ([57|) . with 1 ^ M ^ A^, onto the sub- 
space /r7vJ> and using the refined bound dr(/i)| t-<jv — Ai(//Ar) < C(EI + C) 
(see Remark E]) , we get 

frNHNfrN > fM 



-Xl{HN)gMfi 




C\ 



TM 



frN{^ + C)fr 



N 



(60) 



Here we have used the fact that fufrN = fu when M < rN/2. Finally, we 
can use Lemma [2^ with M = rN, then estimate frNHNfrN by ([60]) . and 
employ the inequality glj < 2J\f+/M < C(m + C)/M and ([#]). We have 



H^ 



N 



> I 



M 



1 





C\ — 



when 1 ^ M ^ A^. Consequently, 

\l{Hn) > Al(IH) - Cl 
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In the latter bound, by choosing M = N^^^, we get 

Xl{Hn) > Al(IH) - C{eR{N) + N-'/^). 
Spectral gap. By (i) one has for every L = 1, 2, 

limmf {Xl{Hn) - Xi{Hn)) = Al(IH) - Ai(]H) 

Taking the limit as L — t- oo we obtain the spectral gap 

liminf (inf aess(^iv) - Xi{Hn)) > inf (TessW - Ai(IH). 

N^oo 

Step 2. Convergence of lower eigenvectors. First we consider the conver- 
gence of ground states. If ^' at is a ground state of i^Ar, then $Ar := Un'^n 
is a ground state of '■= Un{Hn — Neu)U^ on T^'^ . From the proof in 
Step 1, if max{NeR{N),VN} < M < iV, then we have WgM^NW^ and 

lim ifM^NMi^N) = Ai(]HI) = ($(i),EI^.(i)\ 

where ^^^^ is the unique ground state of H on 

To prove <I>Ar — ;> ^>(^\ it suffices to show that /m^n ^''^^ ■ Let us write 

fM^N = aN + bN 
where qn G Span({<I)(^)}) and bN-L^^^^- Then 

{fM^N,MfM^N) = {aN,MaN) + {bN,mN)>Xi{m)\\aN\\^ + X2iM)\\bN\\ 
= ||/M«'iv||'Ai(M) + (A2(M) - Ai(]H))||6^||2. 

Since {fM^NMi^Nh^ - WfM^NW^Xiim ^ and A2(IH) > Ai(]HI), we 
conclude that 5iv— >-OasA^— )-oo. Therefore, /m^n and hence 

$7V ^ ^^^K as N — )• oo. 

Remark on the convergence in quadratic domain. We can show that if (A3s) 
holds true, then we have the strong convergence $Ar — )• ^^^^ in the norm in- 
duced by the quadratic form of IH on J"+ , namely ($Ar,]HI$iv) — ^ (<I'*^^'', IHI$(-'^)) 
In fact, from the above proof we already had 

(/M$iV,IH/M$;v) ^ (^^'\IH«>«) and (<7Af ^7V, ^iv5M$iv) ^ 0. 

On the other hand, from (A3s) and ([IS]), we get Hn > co(]HI + 7V+) + g{N) 
on J-^^ , where cq > and ^(A^) — )• as A — )• oo. By choosing M such that 
max{g{N),N£R{N),y/N} < M < A, we obtain 

guHNgM > co(IH + M) + g{N) > cqM. 

It implies that (gM^N^'^QM^N) — ^ 0. By using the localization for EI in 
Proposition [23] we can conclude that ($Ar,IH$Ar) IH^>(i)) . 

The convergence of excited states of Hn can be proved by using the above 
argument and the following abstract result. 

Lemma 25 (Convergence of approximate eigenvectors). Assume that A is 
a self-adjoint operator, which is bounded from below, on a (separable) Hilbert 
space, with the min-max values Ai(^) < ... < Xl{A) < mfacss{A). If the 
normalized vectors {a^n^}^>i~^ satisfy, for all i,j G {1,2, 

lim = 6ij and lim {xI^\AxI^^) = Xj{A), 
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then there exists a subsequence {xi['f}}k>i which converges, in the norm in- 
duced by the quadratic form of A, to an eigenvector of A with the eigenvalue 
XUA). 

The proof of Lemma[52]is elementary, using the same argument of proving 
the convergence of ground states, and an induction process. The proof of 
Theorem [2] is finished. □ 

7.2. Proof of Theorem El 

Step 1. Convergence of the free energy. We need to show that 
lim {Fp{N) - Nen) = -/^-^og Trj-^ [e-^- 

N—^oo 

We can rewrite 

Ffs{N) - Nen = inf {Tr[^^r] - p-^S{r)} 
r>o,Tr <iv(r)=i 

where Hn := Un{Hn — Ne}i)U'^, and 

-/3"ilogTrj-+[e-H = Trpf] - P~^S{T) 
where T := Z^^e'^^ with Z = Tr [e"'^'^] . 

Upper bound. Let us write T = Xli^i |^^*^)(*^^*'' | where {<I>(*)}j^^ is an 
orthonormal family in and ti > t2 ^ ••• ^ 0, = 1. Then 

oo 

TV^+ [Mr] - r^S{T) = + r\^ogti) . (61) 

i=l 

Fix L £ N. By using Lemma [23] and the fact that EI is bounded from 
below, we can find for every M > 1 a family of normalized states C 
such that limj\//--5>oo(^7i*^> ^j^/ ) = ^ij and 

limsup(M) „ < (]H)^„ (62) 

M-s>oo M 

for all i,j E {1,2, Denote 6l = Yld=i'ti and 



1=1 



Then it is easy to see that Tl^m > and Tr[rL^M] = 1- Moreover, because 
limM^oo(^A],^M) = ^ij we get 



Jim^5(r„,) = -Ej^>°^(Jt) 



L 

(63) 



Choosing M = N^^^ and applying Proposition [15] we obtain 



Tr <N [HnTl,m] - Tr^+ [HTl.m] ^0 (64) 
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when iV oo. From ([63]) and ((Ml), we find that 

hmsup {F,3{N) - Nen) < hmsup fxr <iv [^^yr^.A/] - r^S{TL,M) 



N^oo 



< 



1=1 



Finally, taking L — )• cxd and noting that 9l ^ 1, we obtain from ()6ip that 



limsup(F/3(A^) - A^ch) < -/3"Mog Trj-_^ 

Af-s>oo 



. Us- 



Lower bound. Let us denote Fat := UnT/^^nU^ = e l^^i^ jTi 

ing Tt[H^T ]^] — P~^S{Tn) < and the stability 

TtIRnTn] - /3o"'S(Fjv) > Fp,{N) > -CN 

with 13^^ > /3-\ we obtain T:i[HnTn] < CN. Therefore, by using Lemma 
the simple bound EI < C{Hn + CN)., and Proposition [T5] we find that 



CN 

'W 



IM 



CN 



with F 



<M 
N 



Im^nIm and F 



>M 
N 



9m^n9m- On the other hand, using 



/m + 9ai ~ ^ ^^'^ Brown-Kosaki inequality [lOJ , we have 

s{T^)<s{r^''') + s{r>^''). 

If we choose M = N^^^ , then the above estimates imply that 
F^(A^)-iVeH = (TY[i?^F^] - /3-i5(Fjv)' 

> (l-CiV-^/^)Tr[HF|*^] -/3-i5(F^ , 

+ Tv[HNr^^^] - r'^s(r>^'^) - CN-^i^. 

By using 




and 



we can conclude that 

Fp{N) - Nen 



rp -p>M' 

log(TY F^*^) + ^ ^-'^ ^>^'^ 



TrF 



N 
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Assumption (A3') implies that TrF^^^ — )• and TrF^^^ — t- 1 as — >• oo. 
Therefore, 



liminf (F«(A) - Ach) > hminf ( -/3"^ logTr e 

= -r'logTre-^=. 



Here in the later equality we have employed the Dominated Convergence 
Theorem using Tre"^"''^ < oo for /3o < /3 in (A4). 



Step 2. Convergence of Gihhs states. From the above proof, we have 



where 



lim TrpFjv] - ^ S{Tn) = Tr[MF] - r'SiV) 

A-i-oo 



We have proved that Tv^^pNe-^^^^-^^^'^U*^] Tr^+ e~l^^ as A ^ oo, 
and we will now show that F^v — )• F in the Hilbert-Schmidt norm. It is 
well known that a sequence of non-negative operators Ajsf with Tr{A]\f) = 1, 
which converges weakly-* to an operator A, converges in the trace norm if 
and only if Tr(^) = 1 (see, e.g., [TT], [Ml Cor. 1] and pi Add. H]). Using 
this fact, we will get the result. Moreover, by using Tr[(FAr-FAr)2] ^ due 
to the condensation (A3'), it remains to show that Tr[(FAr - F)^] 0. The 
latter convergence follows from the equality 

TrpFjv] - r^SiVN) - (TrpF] - r^S(r)) = T> [Fjv(logr;v - logF) 

and the following entropy estimate, which is inspired from \27\ Theorem 1]. 
The proof of Theorem [3] is finished. □ 

Lemma 26 (Relative entropy inequality). If A and B are two trace class 
operators on a Hilbert space and 0<j4<1, < B < 1 — e for some e > 0, 
then there exists > such that 

TT[A{logA-logB)] > CeTT[{A- B f] + TiiA- B). 

Proof. A straightforward computation shows that for all x,y £ [0, 1], 

y - 1 - log y 



x{logx- logy) - (x-y) > g{y){x - y) with g{y) :-- 



Since the function y i— )• g{y) is decreasing, we have g{y) > g{l — e) > for 
all ye [0,1 — e]. Therefore, By Klein's inequality [62l P- 330], one has 

Tr [A(log A - log B)] > g{l - e) Tr[{A - Bf] + Tr(yl - B). 



The proof of Lemma [26] is finished. 



□ 
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Remark on stability and condensation. We prove the remark that if (A3s) 
holds true for some ^ (0, 1) and Trjr^ |^g-(i-eo)^oHj ^ f^^. gome /3o, 
then (A3') and (A4) hold true. 

In fact, using M < dr{h + C) + C and Trj-^ [e-^^-^o)*"^] < oo we obtain 
that Trjr^ |^g-(i-eo)/3odr(/i.+C)j ^ ^^^^ ^^^^ f^j. g^gj-y positive operator A 
in we have Trjr^ e-'^^i'^) < oo if and only if Tr^^^ e~^ < oo. Therefore, 
we can conclude that Tr e^^^^^o)^"'* < oo and Trjr_^ e"(^"^o)^o'^^('') < oo. 

Prom the latter bound and the inequality Hj^ > (1 — cq) J2iLi ^« + '^i^) 
(A3s), the stability follows: 



> 
> 



+ 0{N) 



-(l-eo)ft)dr(/i) 



Ffs,{N)-Nen = - 

-/3o-ilogTr_^. 

-/3o-ilogTr^_ ^ 

Next, using Hn > coA/'+ + o{N) on and the stability 

Tr[i?jvr/3,7V] - /3o"'5(r/3,Ar) > F^o(^) " > o{N) 

with r^,Ar := U]\jTi3^]\jU^, we find that 

> Fp{N)-NeH = TT[HNfp,N]-p-'S{fp^M) 



+ o{N) > -C + o{N). 



> 



/5 
/3 



/3 



coTr[AA+r/3,7v] +o(A^). 



Therefore, TrfAA+T^j,^] = Tr[A/'+r/3,7v] = o{N). 



Appendix A. Bogoliubov Hamiltonian 

In this appendix we prove some general properties of quadratic Hamilto- 
nians in bosonic Fock spaces, which includes the results in Theorem [TJ Most 
results of this section are well known [8] . In some parts of the discussion we 
shall refer to [56\ [H] for more details. 

Let 7i be an arbitrary (separable) one-particle Hilbert space. We can 
identify T-L with its dual space 7i* by the conjugate linear unitary operator 
J : Ti ^ H* defined as {Jf){g) = (/, g)n for all f,g £ H. If ^ is a subspace 
of L^(r2), then J is simply the complex conjugate. 

For any vector f £ Ti, we may define the annihilation operator a{f) and 
the creation operator a*{f) on the Fock space J^{Ti) = 0^=o 'S)sym^ ^ 
([37|) and (j38|) in Section HI For any state $ in the Fock space J^(7i), the 
one-body density matrices (7<i>,a#) are operators on V. defined by 

{f,7<s>9) = {a*{g)a{f))^, {f,a^Jg) = {a{g)a{f))^ ior all f, g e n. 

We will be particularly interested in states ^>'s with finite particle expec- 
tation, namely Tr7$ < oo. In this case, (j^^a^) belongs to the set Q, which 
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contains all pairs of operators (7, a) on % such that 7 is trace class, a is 
Hilbert-Schmidt, = a and 



a* 1 + J7 J* 

Let us denote by Qq the set of all pairs (7, a) £ G which satisfy 

aa* = 7(1 + J7J*) and 'ya = aJjJ*. (65) 

The significance of Qo is that for any (7, a) G Qo, there exists a unique 
quasi-free pure state ^> € J^(Ti) such that 

(7, a) = (7$,a$). 

Any element in Q is also associated with a unique state, but the latter is a 
mixed state. 

The one-body density matrices offer a simple way to define quadratic 
Hamiltonians. More precisely, let H he a self-adjoint operator on Ti and let 
K he a Hilbert-Schmidt operator on Ti such that K = and such that 
the following inequality holds true 

A:=(^ ^* JHJ* ) ^ ^ > ™ n®n\ (66) 

We shall consider the quadratic Hamiltonian H on J-'{T-L) defined by 

(M)^ := g(74,,a$) = Tr[Hj^] + 3?Tr[E:a$] (67) 

for every state $ living in the truncated Fock spaces and in the domain of 
H: 

M n 

U 0(2)^(^)- (68) 

M>0 n=0 sym 

It can be verified that the so-defined operator EI is exactly the Bogoliubov 
Hamiltonian given in ()12p . 

The main properties of the quadratic Hamiltonian H are given in the 
following 



Theorem 27 (Bogoliubov Hamiltonian). Let H he defined by ( [<?7| ) and as- 
sume that \66\) holds true. Then we have the following statements. 

(i) (Form domain). We have the quadratic-form inequalities 

C-^dT{H) -C <M< dT{H + C) + C (69) 

As a consequence, the form domain of the Friedrichs extension ofM (still 
denoted by M) is the same as that of dT(H) on Tj^. 

(ii) (Variational principle). For any (7,0:) E Q, we can find (7', a') S Go 
such that 5(7', a') < 5(7,0) and the inequality is strict expect when (7,0) G 
Gq. As a consequence, the ground state energy ofM is 

inf (t(BI) = inf 5(7,0)= inf 5(7,0). 

{■y,a)eg ij,a)ego 

{Hi) (Ground state and ground state energy). The Hamiltonian H has a 
unique ground state in F{H), which is a pure quasi-free state. Moreover, 
we always have inf o"(H[) < except when K = in which case inf a{M) = 
with the vacuum being the corresponding ground state. 
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(iv) (Spectrum). We have aessO^) = cr{M.) + acssiH) and 



ct(H) = inf(T(EI) + < '^riiXi \ Aj G a{SA),{ni} C N and < oo 



i>l 



where 




{v) (Lower spectrum) . // JKJ = K, JHJ* = H and H has infinitely many 
eigenvalues below its essential spectrum, then H also has infinitely many 
eigenvalues below its essential spectrum. 

On the other hand, if JKJ = K > 0, JHJ* = H and H — K has 
only finitely many eigenvalues below its essential spectrum, then EI also has 
finitely many eigenvalues below its essential spectrum. 



Remark 10. Theorem [1\ follows from Theorem 27, with H = h + Ki and 
K = K2. In particular, 06]) follows from / fgp|) because C~^{h + 1) < h + 
Ki < C{h + 1). Moreover, by following the proof of 169\) we also have M > 
dT{h - C) - C; and if Ki = K2 > 0, then M > dT{h - e) - for every 
e > 0. 



Proof. 1. The variational principle is well-known (see, e.g., [441 Theorem 
1.7 p. 101] for a proof). If -fC = 0, then from and the inequality H > r] > 0, 
we see that inf (t(EI) = and the vacuum is the unique ground state of HI 
(which corresponds to 7$ = and 0$ = 0). On the other hand, if 7^ 0, 
then by taking a normalized vector v € T-L such that {Jv, Kv) > and then 
choosing the trial operators 



7A := X\v){v\ ,ax = -^/X(TTXj\v){Jv\ 



with A > small, we can see that inf (t(IHI) < 0. 

Now we show that IHI is bounded from below. If (7,0) G Go, then we can 
find an orthonormal family {n„}„>i for Ti such that 



"f{x,y) = ^ A„n„(x)u„(y), a{x, 

n )Un{x)Un{y) 

n>l n>l 

where A„ > and Yl'^=i An = Tr7 < cxd. Thus we can rewrite 

g(7, a) = ^ iXn {Un, HUn) + \/A„(1 + XnW (u^r, KUr, 



n=l 
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The inequality ([66|) implies that {un,Hun) > \ {un,Kur, 
Therefore, 



+ r] for all n > 1. 



n=l 

oo 



> r7Tr7-^ 



^/XJT+Xn) ) I {Un, KUr, 
An I (^i^, KUn) I 



n=l 



> 



> 



Tr 7 




(70) 



r7Tr7 - V ±r7||j\ ||hs, 

where ||-fC||HS is the Hilbert-Schmidt norm of K. Thus EI > —C. 

2. Next, we show that C-^dT{H)-C < EI < dT{H+C)+C. In fact, from 
the above result, we see that the quadratic Hamiltonian with A replaced by 

H-r]/2 K 
K* JHJ*-r]/2 

is also bounded from below. Therefore, 

C. 



>7]/2 



= {ri/2) Tr[7$] + Tt[{H - 7?/2)7$] +^Tr[Ka^] > {ri/2) Tr[7$] 

Similarly, for a constant Co > large enough, one has 

TV[C7o7$] + RTviKa^] > -C and Tr[C7o7$] - K Tr [i^a$] > -C. 

These estimates yield the desired inequalities. 

3. Now we show that EI has a ground state. Let a sequence {(7^, an)}$?Li ^1 
Qo such that 

lim q{jn,an) =infcr(IHI). 



The inequality (f70]l implies that Tvyn is bounded, and hence Tr(a„a*) is 
also bounded. Thus up to a subsequence, we may assume that there exists 
(70 ) 0(o) £ Q such that ^ oq and 7„ ^ 70 weakly in the Hilbert-Schmidt 
norm. Consequently, lim„_!.oo Tr[/Ca„] = Tr[Erao] and liminf„_>.oo Tr[ii'7„] > 
Tr[i77o] by Fatou's lemma since H > 0. Thus (7o,Q!o) is a minimizer of 
q{'y,a) on Q. Due to (i), this minimizer (70,00) belongs to Qo- 

4. To understand the structure of one-body densities matrices and the 
spectrum of the quadratic Hamiltonian, let us introduce Bogoliubov transfor- 
mations. A Bogoliubov transformation V is a linear bounded isomorphism 
on such that (W* — 1) is trace class (Stinespring condition) and 



V 




J 



J* 





J* 
J 



V, VSV* =S :-- 





-1 



Since (70, ao) £ Go, there exists a Bogoliubov transformation Vo on HQH* 
which diagonalizes the one-body density matrices (70,0:0), namely 



Vo 



To 

On 



Oo 

1 + J 10 J* 




1 



(see, e.g., [561 E]). Then by employing the fact that (70,00) is a mini- 
mizer for 5(7, a) on G, we can show (see [44^ Theorem 1.7 p. 101] for a 
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detailed proof) that the Bogohubov transformation V := SVqS diagonahzes 
A, namely 

= ( jIj- ) 

for some operator onV.. From ([66|) and (f7T]l . we deduce that ^ > > 0. 

Using the Bogoliubov transformation V, we can find a unitary transfor- 
mation Uy on the Fock space F{T-L), called the Bogoliubov unitary (see e.g. 
[561111]), such that 

oo N 

UvHUt, - inf a(H) = dr(e) := J^^i- 

7V=0 i=l 

This implies that EI has a unique ground state and ct(EI) = inf o"(IH) + dr(^). 

5. Next, from (f7T]l and the assumption VSV* = S, it follows that 

J^J* ) ~ ^ " SVAV* - A = SVS{SA - X)V*. 
Thus (t(^) = a{SA) because V is an isomorphism on T-lQH*. Consequently, 
cj(EI) = inf cr(EI) + I ^n^Ai | Xi G a{SA),{ni} C N and < ooj. 

i>l 

6. Now we prove that o"ess(0 = cress{H). From (f7T|) . it follows that 
^iC) = cj(V^V*). On the other hand, since K is a compact operator on Ti, 
one has cJess(-4.) = (7css{H). On the other hand, by using the identity 

V^V* - A = V{A - A)V* + A(VV* - 1) 

and the fact that (VV* — 1) is compact (indeed it is trace class), we obtain 
cJess(V^V*) = C7ess(^)- Thus o-ess(0 = cTcssiH). Consequently, from cr(EI) = 
inf o-(]H) + cT(dr(0) we obtain 

aess(IH) = inf a(M) + aess(dr(0) = ^(H) + aessiH). 

7. From now on, we assume that JKJ = K and JHJ* = H. In this 
case, (j66p implies that H + K > r] > and H — K > r] > 0. Before going 
further, let us give a simple characterization of the spectrum of H. We only 
deal with eigenvalues for simplicity. 

Assume that t > is an eigenvalue of Then from the above proof, we 
see that t is an eigenvalue of SA. Using the equation 

we find that 

({H + K)x = ty, 
\{H- K)y = tx, 

where x = u + v and y = u — v. Thus 
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Note that x^Oifu(Bv^O. Therefore, is an eigenvalue of 

In fact, using Weyl sequences and arguing similarly, we also obtain 

t e a{^) = a{SA) ^Oe a{Xt). (72) 

8. We show that the number of eigenvalues (with multiplicity) of ^ below 
fi := inf (Tess(0 is equal to the number of negative eigenvalues of X^. 

Note that the mapping 1 1— ?■ At is strictly decreasing. As a consequence, for 
every j = 1,2,..., the min-max value Aj(Aft) is a continuous and decreasing 
function on t > 0. More precisely, if < ti < t2, then Aj(AftJ > Xj{Xt2) 
and the inequality is strict if Aj(A'tj) is an eigenvalue of A'j^. Moreover, 

inf a{Xo) = mi a{H + K) > r] > 

and, for every t G (0, /u), 

inf (Tess (^i) > inf W = 0. 

Therefore, by using ()72p . we obtain a one-to-one correspondence 

(7(e) n (-oo,/i) o (7(^-^)0 (-00,0). (73) 

9. From the inequality fi^{H - K)'^ + H - K > 2fi we get 

Moreover, note that inf cJessC'^/x) = inf (T(.ss{H — h) = 0. Therefore, the num- 
ber of negative eigenvalues (with multiplicity) of is not less than the 
number of negative eigenvalues (with multiplicity) of {H — /x) . 

In particular, if (H — /j,) has infinitely many negative eigenvalues, then 
has infinitely many negative eigenvalues. By ()73p . we see that ^ has 
infinitely many eigenvalues below its essential spectrum. Consequently, H 
has infinitely many eigenvalues below its essential spectrum. 

10. Now we assume furthermore that K > 0. Then we get the inequality 

Moreover, note that 

inf aessiX^) = inf a^ss (^H - K - ^ = 0. 

Thus, if H — K has finitely many eigenvalues below ^u, then Xt also has 
finitely many negative eigenvalues. By (j73p . ^ has finitely many eigenvalues 
below its essential spectrum. Consequently, IH has the same property. □ 

Appendix B. Localization of band operators on 7"+ 
In this appendix we prove the localization in Proposition [22l 
Proof of Proposition Using the IMS-identity 

A = fkiAfM + ghiAgM + ^ ([/m, [fM,A\] + [qm, [qm^A]]) 
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we can write, with $ = e°^o^i ^ ®i^o-^+ = -^+' 
{<5,A^) = {fM<^,AfM^) + {gM^,AgM<^) 

+ E [(^/ (i) - fM (j)? + {gu (i) - gM {^^,A^J). 

l<\i-j\<(T 

Since / and g are smooth, we have 

ifM (0 - fM {3)f + {gu ii) - gu {j)f < (ll/'lloo + ll^'lloo) ^^J^. 

Moreover, using the assumption that ^ > we get 

l<|i-il<o" l<l*-j|<o- 

i 

Therefore, we obtain the operator inequahty ()54p . 
Finally, let us show that if 

S := sup{||ffAf^|P : $ G y, ll^ll = 1} < (dimy)-\ 

then dim(/M^) = dimF. In fact, assume that dimF = L and let an 
orthonormal basis for Y. For all {ai}^^^ G C^\{0}, we have 

L 2 

i=l l<i<j<L 



i=l 
L 

= Y\(^i\'^('^-\\9M'^if)-'^ J2 \'^iaj\-\{gM^i,gM^j)\ 

i=l l<i<i<L 

L L 

i=l i<i<j<L i=l 

Therefore, the subset {fM^j}j=i C /mY is linearly independent, which 
implies that dim(/Afy) > L = dimy. On the other hand, dim(/7v/l') < 
dim!" because /m is linear. Thus dim(/My) = dimY. □ 

Appendix C. Logarithmic Lieb-Oxford inequality 

We provide a proof of Proposition [TTJ We follow ideas from [33| [371 [3H]. 

Proof of Proposition \11\ . Let us write for simplicity p := p^,. Note that for 
any fi^i ^ which is radially symmetric about Xi and such that / fixi = 1, 
we have 

w ^ ^ 

\ i=l i=l / 

1 ^ N N 

= l^D{p,p)+ E D{px^,px.) + -YD{px,,^^x^ -Yl^^P'^^''^^- 



2 

l<i<j<N i=l i=l 
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We have used here Proposition[8]together with the fact that J^2{p—YliLi f^x,) 
0. By Newton's theorem, {fixi * w){x) < wi^x) a.e., and therefore 

D{lJ,x,,tJ'Xj) < w{Xi - Xj) 

for any i ^ j. We arrive at the following Onsager-type estimate |45] 

1 ^ 1 ^ 

w{xi- Xj) > --D{p,p) + YD{p,fix,) - - ^Di^ix^^fJ^xJ- (74) 

l<i<j<N i=l 1=1 

Writing D{p, fi^i) = {p*'w){xi) + D{p, fi^i — SxJ and taking the expectation 
value against ^, we get 

( Y] w{xi-Xj)) >\d{p,p)+( p{x) D{p,px - Sx)dx 

- o / p{x)D{lix,Px)<^X. 

Now we choose 



p,{y) = R{x)-^ Jy—^\ and R{x)- '^^^^ 



R{x) ) y p{x) ' 

where p = 7r~-'^x(|x| < 1) is the (normalized) characteristic function of the 
unit ball. We will choose the function A(x) at the very end. 
We start by computing 

f f 1 p(x) 

D{nx,Px) = - / p{y) p{z) log {R{x)\y - z\) dy dz = -log— -+D{p,n) 

which gives 

-\ \ p{x)D{px,Px)dx = -] I plogp-l-D{p,i^i) I p + l I plogA. 
Hence we have proved that 

( V w{xi-Xj)) >]-D{p,p)-- I plogp-]-D{p,p) I p 

\l<t^<N 2 Ah2 2 J^2 

+ - / plogA+ / p{x) D{p,p.x - ^x)dx 
4 Jm? Jr2 

and it remains to bound the last term. 

Let vr = {7rR)~'^x{\x\ < R) be the (normalized) characteristic function 
of the disk of radius R. We have 

1 / IxP IxP \ 
{uR*w){x) = w{x) - 2 ( ~^°S"r2" + - M ^(1^1 - ^) 

and thus 

{uR * w){x) - w{x) > log ^J^^ x{\x\ < R). 
By scaling we find that 

D{p, px - Sx) = D{p{- - x),vfi(x) - ^o) > / p{y - x) log ( 4rT ) 

J\y\<R(x) \R{X)J 
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and therefore 

/ p{x) D{p,fi^ - 6x)dx > / dxp{x) dyp{y-x)log 
JIR2 JiR2 J\y\<R{x) 



R{x) 



The followmg is similar to [331 Lemma 2]. 
Lemma 28. Let f G Li(IR^M+) and let 

M/:=sup^/ f{y)dy. (75) 
r Trr J\y\<r 

Then we have for any R > 



J ^^dy/(y) log I^^MfR". 



\y\ 



(76) 



Proof of Lemma [2R Let /(r) := f^^ f(r cos (6), r sm{d)) dO be the spherical 
average of /. Then we have 

- dy f{y) log j ^ ~ / '^^^^^ (^) '^^ 

< vrM/ / rdr = -MjR^. 
Jo 2 

In the previous estimate we have first integrated by parts, and then used 
J^sf{s)ds = J^J<7rr^Mf. □ 



Using the previous estimate (j76p with x G M fixed and recalling our 
choice R{x) = ^JX{x)/ p{x), we get 

Aypiy - X) log ( Jj^) > -| M,(x) R(xf = M,(^-) 



/ 

J\v\ 



\y\<R{x) 
where 



Mp{x) := sup [ p{y - x) dy 
r 7rr^ J\y\<r 

is the Hardy-Littlewood maximal function of p. This gives the estimate 

/ dxp{x) f dy p{y - x) log ( ^1^] > f 

Jr2 J\v\<r(x) \^\X)J ^ Jm 



and our final bound is 



V w{xi-Xj)) >\D{p,p)-\f plogp-^D{p,p) [ p 



Here the optimal choice is A = p/{2-KMp) but the term Jj^aplogMp is 
then not easy to control. At least we can take A = (p/ J pY and we find an 
error term of the form 



3/ "'"^-f-^i [t-^' 

4 7r2 } P 2 J^z VJ P 



p 
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which we can also write as 



p]\--J /log/-^ / fMf), where/ 



Now we use that ef log / < /(I + /^) and that / < Mj to get an error of 
the form 



Using now 



2 ./to2 ■' t 



we end up with a total error of the form 

C 



□ 
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